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ABSTRACT  .}/ 

V<  ^ ' t\ 

Experiments  on  the  propagation  of  narrow,  high  amplitude  pulses 
along  a plasma  column  have  been  conducted.  The  evolution  of  the  pulse 
in  space-time  has  been  observed  as  a function  of  initial  pulse  strength. 

This  technique  gives  insight  into  the  competitive  effects  due  to  dispersion 
and  nonlinearities  . 

For  a plasma,  the  positive  column  of  a hot  cathode,  D.  C.  discharge 

in  Argon  is  used.  The  pressure  varied  from  1 to  10m  iTorr  Jig,  and  the 

9 i 10  A3 

averaged  electron  densities  are  of  the  order  of  10  f-  10  ctrr  . The  pulses 
for  excitation  are  generated  by  a home  made  device  capable  of  producing 
subnanosecond  pulses  with  3200  volt3  peak  amplitude,  into  a 50C  line,  at  a 
repetitive  rate  of  approximately  lOOp.p.s.  ^ ^ 

At  each  operating  point,  i.e.  , for  a given  pressure  and  discharge 
current,  recordings  were  obtained  for  different  values  of  excitation  pulse 
amplitude  and  background  static  magnetic  field. 

For  zero  static  magnetic  field,  two  types  of  recordings  were  taken. 

For  the  first  type,  the  position  of  the  receiver  was  fixed  and  oscillograms 
in  time  were  recorded.  In  the  second  type,  the  sampling  time  is  kept  fixed 
and  the  receiving  probe  is  moved  along  the  column.  From  the  oscillograms, 
we  observe  that  the  dispersion  is  negative  and  that  tlie  group  velocity  increases 
with  plasma  density.  As  the  amplitude  of  the  input  pulse  is  increased,  non- 
linear effects  start  to  play  a role.  Under  experimental  conditions,  these 
effects  did  not  lead  to  quantitatively  different  behavior  than  linear  phenomena. 
From  numerically  taken  Fourier  transforms  of  the  observed  propagating 
wavepackets,  the  transfer  of  energy  between  different  frequency  components 
forming  the  disturbance  have  been  observed.  The  region  of  interaction  in 
the  frequency  domain  is  seen  to  be  a function  of  distance  along  the  column 
and  of  wave  velocity. 

For  finite  static  magnetic  field,  recordings  of  the  first  type  were 
obtained  as  a function  of  pulse  amplitude.  Again,  the  dispersion  is  negative 
and  the  group  velocity  increases  with  both  and  u>c-  Contrary  to  the  Bq=  0 
case,  for  large  pulse  amplitude,  qualitatively  different  wave  phenomena  were 
observed:  in  particular  the  nonlinear  waves  show  modulational  splitting. 

This  phenomenon  was  observed  both  in  time  and  in  the  frequency  domain. 
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I.  Introduction 

The  general  problem,  we  are  concerned  with  is  that  of  the  inter- 
action of  a plasma  medium  with  an  electromagnetic  field.  At  linear 
level  of  excitation,  the  various  problems  of  wave  propagation  have  been 
explored  at  length  and  a fair  amount  of  understanding  has  been  reached. 

At  nonlinear  levels  of  excitation,  due  to  the  lack  of  a general  method 
of  attack,  the  analysis  made  to  date  are  mostly  ad  hoc  and  under  various 
forms  of  idealization  and  approximation.  The  purpose  of  this  investiga- 
tion is  to  study  some  of  these  interactions  both  theoretically  and  ex- 
perimentally. The  results  have  yielded  new  data  and  insight  into  the 
problem  of  the  competition  between  dispersion  and  non-linearity  in  a num- 
ber of  plasma  configurations.  The  comparison  between  theoretical 
result  and  experimental  data  is  generally  favorable. 

Theoretical  studies  in  the  linear  and  weakly  nonlinear  regimes 
may  be  grouped  into  two  categories,  namely,  analysis  in  'unbounded 
and  bounded  plasmas.  Experimental  investigations  on  the  other  hand 
must  be  performed  with  plasmas  of  bounded  geometry.  There  are  ex- 
periments, however,  where  the  pertinent  dimensions  are  such  that  to 
a good  approximation,  the  waves  under  consideration  may  be  thought 
as  propagating  in  an  infinite  medium  and  the  results  behave  as  predicted  by 
unbounded  plasma  theory.  ^ Studies  in  unbounded  plasmas  have  been 
many,  and  several  text  books  have  been  written  on  this  subject.  Com- 
prehensive expositions  have  been  done  by  Ginzburg^  and  Stix^  for 

(4) 

the  linear  regime;  and  for  the  weakly  nonlinear  regime  by  Tsytovich, 

Davids  on  and  F.  Einaudi  et.  al. 

The  bounded  geometry  considered  in  this  investigation  are: 
a circular  plasma  column  imbeded  in  an  infinite  dielectric  (free  space) 
and  a circular  magneto  plasma  column  surrounded  by  a conducting 
waveguide  of  the  same  dimension  as  the  plasma  column.  These  struc- 
tures are  known  to  support  numerous  types  of  waves.  The  work  in 
this  study  is  concerned  with  the  propagation  of  the  slow  mode  whose 
characteristics  do  not  depend  on  ion  motion.  These  waves  are  some- 
times called  "space  charge"  or  "electron  plasma"  waves.  They  differ 
from  the  "true"  electron  plasma  waves  in  an  infinite  plasma  since  the 
boundary  conditions  posed  by  the  system  have  altered  the  wave 
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structure.  Practical  interest  in  these  modes  are  the  possibilities  of 
aoolications  in  laboratory  plasma  diagnostics'1’  \ high  gain  micro- 
wave  devices  , and  plasma  heating 

The  basic  experimental  results  of  linear  slow  wave  propagation 
have  been  well  documented  ^ and  good  agreements  with  theoret- 


ical analysis  have  been  obtained 


15, 16) 


However,  the  situation  is 


different  for  the  weakly  nonlinear  regime.  Very  few  theoretical  analy- 

(17  18  19  ?0 ) 

sis  of  the  boundary  value  problem  have  been  done  ’ 1 ’ The 

emphasis  in  most  of  these  studies  has  been  on  the  nonlinear  coupling 

(21) 

of  waves,  in  particular  three  wave  interaction  . This  problem 

(22) 

takes  a variety  of  forms.  Perulli,  et.  al.  studied  the  decay  of  a 
slow  wave  into  two  other  slow  waves  with  different  azimuthal  variation. 
Larsen  used  the  method  of  the  averaged  Lagrangian  to  analyze 
the  interaction  of  three  slow  waves.  Kuhn^-^  formulated  the  same 


problem  from  a coupled  mode  formalism.  Aside  from  experiments 
on  surface  wave  echo  (3  patial)  experimental  work  has  been  confined 
to  the  mixing  of  two  pump  waves  of  different  frequencies  to  produce  a 

/ p L ) 

third  at  the  beat  frequency 

All  of  the  theoretical  and  experimental  efforts  mentioned  thus 

far  have  one  feature  in  common:  they  are  all  concerned  with  steady 

state  situations.  Recently,  transient  techniques  have  been  employed 

in  the  experimental  investigation  of  wave  propagation  in  plasmas. 

Voltage  steps  and  pulses  (base  band  and  RF)  have  been  used  to  excite 

transient  waves.  Schmitt  was  one  of  the  first  experimentalists  to 

use  pulse  excitation  to  study  plasmas.  He  observed  the  dispersion  of 

base  band  pulses  propagating  through  a plasma  filled  coaxial  line. 

( 2.R)  ( 29 ) 

Proni,  et.  al.  and  Tregruis,  et.  al.  used  microwave  pulses  to 

excite  transient  electron  plasma  waves  in  systems  where  the  boundary 

had  no  effect.  The  only  reported  applications  of  baae  band  pulses  to 

study  the  propagation  of  space  charge  waves  were  by  Anicin,  et.  al. 

(31) 

for  the  symetric  mode,  and  by  Demokan,  et.  al.  for  the  dipole 

mode.  Recently  Landt  et.  al.  used  various  types  of  transient  inputs  to 

( 32) 

study  the  linear  properties  of  these  modes  . These  experiments 
were  limited  to  the  linear  regime. 

In  the  nonlinear  regime,  pulse  excitations  and  time  domain 
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observations  have  been  reported  by  Sindoris,  Cheo  and  Grody'  in  a 
Tonks  - Dattner  type  structure  using  a home  developed  device  called  a 
Bouncing  Ball  Generator  (BBG),  which  produced  baseband  pulses  129 
picosecond  wide  at  1.  1 KV  into  50  ohm.  Other  nonlinear  transient  work 
includes  the  observation  by  Ikezi  et.  al.  of  the  propagation  of  elec- 
tron plasma  waves  excited  by  voltage  steps,  and  the  observation  by 
K.  Saeki  et.  al.  of  electron  plasma  wave  shocks.  Manheimer '43* 
has  studied  the  development  of  finite  amplitude  electron  plasma  waves 
in  a bounded  plasma  with  infinite  magnetic  field  using  the  method  of 
expansion  in  terms  of  a complete  set  of  linear  solutions.  He  did  not 
include  the  effect  of  dispersion,  but  predicted  the  steepening  observed 
in  (48). 

The  work  in  this  effort  consists  of  a comprehensive  investigation 
of  linear  and  nonlinear  phenomena  affecting  the  propagation  of  very  short 
base  band  pulses  along  a glow  discharge  magneto  plasma  column.  The 
advantages  of  this  approach  over  nonlinear  steady  state  studies  are:  1) 
by  keeping  the  duty  cycle  low,  heating  of  the  background  electrons  (a 
major  problem  in  CW  experiments)  need  not  be  considered;  2)  if  the 
transient  response  is  short  compared  to  ionization  times,  background 
electron  density  changes  will  not  occur;  3)  the  ability  to  observe  short- 
lived phenomenon  is  greatly  enhanced  because  of  the  expanded  time  re- 
solution. The  major  advantage  of  transient  studies  both  large  and  small 
amplitude  excitations  is  that  we  can  observe  the  development  in  time  or 
in  space  of  an  initial  disturbance  produced  at  some  point.  Tremendous 
insight  is  gained  into  the  competitive  effects  of  dispersion  and  nonlin- 
earity. 

The  source  that  was  used  to  excite  the  transient  waves  is  an  im- 
proved version  of  the  BBG  mentioned  in  the  experiment  of  Sindoris 
et.  al.  The  new  BBG  produces  a pulse  0.  50  nanosecond  wide  at 

3.  2 KV  into  50  ohm.  This  generator  and  the  techniques  developed 
presented  us  with  the  unique  tools  for  such  comprehensive  studies. 

In  chapter  II  of  this  report  the  experimental  program  is  presented. 
To  analyze  the  experimental  observations  a nonlinear  theory  for  wave 
propagation  along  plasma  column  was  developed.  This  theoretical 
analysis  is  presented  in  chapter  III.  In  Chapter  IV  the  experimental 

results  are  then  analyzed. 

•d> 
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II.  Experimental  Program 

The  basic  experimental  set  up  is  shown  schematically  in  figure 
2.  1,  where  the  plasma  is  the  positive  column  of  a glow  discharge  con- 
fined in  a glasss  tube  175  cm  long,  with  A and  K the  anode  and  cathode 
respectively.  M-M'  represents  a set  of  coils  providing  an  axial  D.  C. 
magnetic  field  up  to  1.  2K  Gauss.  WG  is  a removable  conducting  wall 
surrounding  the  plasma  tube.  L is  a pair  of  parallel  plates  connected 
to  the  pulse  generator  (BBG)  which  establishes  an  impulse  like  elec- 
tric field  in  the  plasma.  The  wave  evolution  in  space-time  as  a func- 
tion of  initial  pulse  strength  is  monitored  by  the  receiving  structure  R 
on  a sampling  scope.  The  work  was  done  in  two  stages.  The  first 
consisted  of  a comprehensive  investigation  of  linear  and  nonlinear 
phenomena  affecting  the  propagation  of  pulses  along  the  positive 
column  of  a glow  discharge.  The  second  stage  involved  the  addition  of 
a uniform  longitudinal  magnetic  field  variable  up  to  1.  2K  Gauss  and  . 
the  surrounding  of  the  plasma  column  by  a conducting  wall.  Subse- 
quently, a description  of  the  experimental  apparatus  and  procedure  is 
given. 

2.  1 Description  of  the  Apparatus 
2.  1.  1 the  plasma 

For  a plasma,  the  positive  column  of  a hot  cathode,  giow  dis- 
charge in  argon  is  used.  Fig.  2.  2 shows  a layout  of  the  discharge  tube 
and  its  associated  vacuum  and  electrical  systems.  The  discharge 
tube  comprises  two  Western  Electric  mercury  vapor  cathodes  con- 
tained in  a round  flask  and  connected  by  means  of  a quick  glass  to  glass 
joints-  to  a 175  cm  long  glass  tube,  the  end  of  which  is  terminated  with 
a hallow  anode.  Because  of  the  quick  coupler,  glass  tube  sections  of 
different  diameter  may  be  used. 

It  was  necessary,  due  to  the  type  (sampling  technique)  and 
amount  of  measurements  that  had  to  be  performed,  to  provide  a clean 
and  stable  discharge.  By  allowing  Argon  gas  to  flow  through  the  dis- 
charge tube  at  a slow  rate;  contamination,  due  mainly  to  ion  bombard- 
ment of  the  cathode,  is  reduced  since  such  impurities  are  constantly 
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removed  from  the  tube.  To  provide  effective  regulation  of  the  gas 
pressure,  the  argon  inflow  is  regulated  by  a flow  meter  valve  and 
needle  valve  and  the  evacuation  rate  by  a high  vacuum  valve.  Adjust- 
ment of  these  valves  provide  the  desired  pressure. 

Once  the  discharge  is  operating,  the  desired  range  of  electron 
densities  is  obtained  by  varying  the  current  through  the  pentode  re- 
gulator circuit.  Stability  of  the  system,  over  all  range  considered, 
was  excellent. 

Basic  parameters  of  the  positive  column,  i.  e.  electron  tempera- 
ture, density,  and  collision  frequency , were  measured.  Electron 

temperature  is  obtained  using  the  well  known  method  of  Langmuir 
( 35) 

probe  . To  measure  average  electron  number  density,  the  cavity 

(9) 

perturbation  method  is  used  . The  cavity  is  mounted  in  a carriage 
which  can  move  along  the  plasma  column.  The  electron  density  varia- 
tion along  the  column  is  found  to  be  less  than  0.  57a.  . During  the  ex- 
periments, the  density  is  constantly  monitored  to  prevent  drifts  from 
the  operating  point,  due  primarily  to  pressure  variations. 

The  cavity  is  also  used  to  measure  the  total  collision  frequency 

of  the  electrons.  By  measuring  the  change  in  Q of  the  microwave 

( 36  ) 

cavity  due  to  the  plasma,  this  parameter  may  be  calculated  . The 

collision  frequency  of  electrons  with  ions  and  with  neutrals  is  also 

(37) 

calculated  using  the  equation 


I - 


discharge  tube,  associated  vacuum  system  and 
circuitry 
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Table  2.  1 

Typical  Discharge  Characteristics 


Radius  of  column 
Neutral  density  (at  1/j) 

Percent  ionization 

Electron  neutral  collision  frequency 

Electron  neutral  collision  time 

Debye  length 

Electron  thermal  speed 

Electron  density 

Electron  cyclotron  frequency 

Electron  plasma  frequency 


. 66  cm 

13  -3 

3.  3 x 101-5  cm 


.01% 

17.2  MHz 


1 00  n sec. 
1 mm 

1 . 5x10®  cm/ 
2 x 10^  cm 
. 632  GHz 


sec . 
-3 


1 . 3 GHz 
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Typical  parameters  obtained  using  the  experimental  techniques 
mentioned  above  are  shown  in  table  2.  1. 

2.  1.  2 The  pulse  generator 

A Bouncing  Ball  Generator  (B3G)  is  the  source  for  the  pulses 

used  in  the  experiment.  This  generator  has  previously  been  used  in 

time  domain  studies  of  EM  precursors  and  in  impulse  stimulated 

( 33) 

emission  from  plasmas  . The  generator,  at  the  time  of  those  ex- 
periments, produced  base  band  pulses  with  peak  voltage  of  1.  1 kilo- 
volts and  a risetime  faster  than  120  picoseconds.  For  this  experi- 
ment, the  B3G  was  modified  and  is  now  capable  of  producing  pulses 
with  peak  voltage  of  3.  2 kilovolts  with  approximately  200  picoseconds 
rise  time.  The  peak  power  into  a 50  ohm  line  is  259  kilowatts.  The 
average  power,  however,  is  4 milliwatts  because  of  the  short  pulse 
duration.  The  display  of  the  BBG  output  pulse  shown  in  Fig.  2.  3 is 
an  X-Y  recorder  plot  of  the  oscillogram  from  the  sampling  scope. 

To  obtain  such  an  oscillogram,  the  pulse,  after  being  attenuated  63  dB 
is  applied  to  the  scope  through  a 60  nanosecond  delay  cable  (RG-9B/u) 
of  4.  4dB  insertion  loss.  The  output  of  the  scope  is  then  used  to  drive 
the  X-Y  recorder.  Also  shown  in  fig.  2.  3 is  the  voltage  spectrum  of 
the  pulse.  This  is  obtained  from  Fourier  analysis  using  a computer  of 
the  time  domain  signal.  The  spectrum  is  extremely  wide,  extending 
from  D.  C.  and  almost  flat  to  1 GHz. 

To  control  the  amplitude  of  the  pulse,  wide  band  attenuators 
(GR  type  874-GL)  are  used  at  the  output  of  the  BBG.  Also  at  the  out- 
put, a pre-trigger  pick  off  has  been  installed  to  obtain  a trigger  signal 
for  the  3cope. 

2.  1.  3 Parallel  plate  structure 

To  couple  the  BBG  pulse  to  the  plasma,  a section  of  parallel 
plate  transmission  line  is  used.  The  parallel  plate  coupling  structure 
consists  of  a wide  band  coaxial  to  parallel  plate  transition  which  opens 
up  to  a length  of  uniform  transmission  line,  then  tapers  down  again  to  a 
coaxial  line.  The  line  is  then  attenuated  by  20  d B and  terminated  with 
a 50  12  load.  The  dimensions  of  the  structure  are  such  as  to  keep  a 
constant  50  ohms  characteristic  impedance  along  the  structure  over  a 
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wide  range  of  -requencies.  i ne  plasma  tube  is  inserted  oerpendicular 
to  the  plane  of  the  plates  through  holes  made  at  the  center  of  the  plates, 
as  shown  in  fig.  2.  4.  The  holes  are  made  to  fit  tightly  around  the 
column. 

The  plate  structure  is  a very  good  wide  band  coupler  with  an  al- 
most flat  passband  extending  from  50  MHz  to  approximately  2 GHz. 
This  feature  of  the  structure  made  it  a suitable  coupler  in  the  pulse 
experiments  because  of  the  wide  bandwidth  of  the  exciting  pulse.  A 
similar  structure  is  also  used  to  couple  the  propagating  wavenacket  to 
the  sampling  scope.  The  receiving  plates  are  mounted  on  a carriage  to 
allow  for  movement  along  the  column. 

2.  1.4  Sampling  scope  and  recorder 

To  observe  the  input  pulse  and  the  propagating  wavepackets,  a 
Tektronix  564  storage  scope  with  sampling  plug-ins  is  used.  To 
facilitate  further  analysis,  the  observed  oscillograms  are  also  re- 
corded using  a Hewlet- Packard  X-Y  recorder  connected  to  the  output 
of  the  storage  scope. 

It  must  oe  remarked  that  the  sampling  technique  for  observing  fast 
transient  time  phenomena  requires  that  each  experiment  be  identical. 
This  in  turn  demands  that  the  input  pulse  be  identical  for  each  experi- 
ment and  that  the  plasma  relaxes  to  its  initial  state  before  the  follow- 
ing pulse  arrives  at  the  launcher.  The  last  condition  may  be  ascer- 
tained by  comparing  the  width  of  the  packet  with  respect  to  the  pulse 
separation.  However,  considering  that  the  BBG  is  an  electromechani- 
cal device,  it  is  remarkable  that  the  pulses  it  produces  are  identical 
within  a few  percent  of  each  other,  to  the  extent  that  the  technique 
was  succesful. 

2.  2 The  Experiment 

2.  2.  1 Zero  axial  magnetic  field 

For  the  investigation,  the  apparatus  described  in  section  2.  1 
was  arranged  as  shewn  in  fig.  2.  5.  The  output  pulse  of  the  BBG,  after 
attenuation  to  the  desired  voltage,  is  coupled  to  the  plasma  through  the 
wide  band  parallel  plate  structure.  A 60  nanosecond  delay  cable  is 
used  between  the  generator  and  the  plates.  For  high  amplitude  pulses, 
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difficulty  is  encountered  with  reflections  from  the  plates  (due  to  a 
small  mismatch),  that  travel  back  to  the  BBG.  At  the  BBC-,  these  are 
once  again  reflected  (the  B3G  is  an  open  transmission  line)  and  re- 
turn to  the  plates  generating  unwanted  signals.  3v  introducing  the  de- 
lay line,  the  reflection  arrives  at  the  plates  120  ns  after  the  original 
pulse,  long  after  the  observations  are  made. 

The  plasma  column  is  inserted  perpendicular  to  the  plane  of  the 
plates  through  the  holes  at  the  center.  The  holes  are  made  to  fit 
tightly  around  the  column.  As  the  pulse  propagates  along  the  parallel 
plates,  the  region  of  plasma  within  the  plates  ( 1 cm  separation!  feels 
uniformly  around  the  column  the  effect  of  the  electric  field  of  the  pulse. 
The  field  is  directed  parallel  to  the  axis  of  the  column,  such  that  elec- 
trons are  accelerated  towards  the  anode. 

With  this  mode  of  coupling,  circularly  symmetric  waves  are 
excited  in  the  plasma  and  they  propagate  along  the  column  towards  the 
anode.  To  prevent  waves  from  propagating  towards  the  cathode,  the 

section  of  olasma  between  cathode  and  launcher  is  surrounded  with  a 

* • 

copper  sheet  and  no  wave  can  propagate  under  those  conditions. 

More  will  be  said  on  this  in  chapter  III.  Since  the  tube  is  of  finite 
length,  wideband  microwave  absorbers  are  used  at  the  end  of  the 
column  to  prevent  possible  reflections  from  the  anode  and  into  the 
receiver. 

The  receiver  is  mounted  on  a carriage  and  can  be  moved  along 
the  whole  length  of  the  column.  These  plates  pick  up  the  longitudinal 
electric  field  associated  with  the  propagating  wave.  The  field  strength 
is  displayed  on  the  sampling  scope  and  also  plotted  on  paper  using  the 
X-Y  recorder.  A trigger  pulse  from  the  BBG  is  used  to  properly 
synchronize  the  time  of  sampling. 

Two  types  of  oscillograms  were  recorded.  First,  for  a fixed 
position  of  the  receiving  plate,  oscillograms  in  time  were  taken.  The 
analog  output  of  the  scope  drives  the  Y axis  of  the  plotter,  while  the 
scope:s  time  base  is  used  to  drive  the  X axis.  In  this  manner,  a replica 
of  the  trace  that  appears  on  the  scope's  screen  is  plotted.  Secondly, 
the  spatial  distribution  of  the  waves,  for  fixed  sampling  time,  were 
obtained.  The  Y axis  of  the  plotter  is  driven  as  in  the  time  measure- 
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Mock  diagram  of  impulse  generation,  transmission  and  reception 
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rr.ents;  the  X axis  however  is  driven  by  a voltage  proportional  to  the 
position  of  the  receiver  .vith  respect  to  the  launcher.  The  arrangement 
is  shown  in  fig.  2.  6.  A constant  voltage  is  applied  across  the  potentio- 
meter. As  the  crank  is  turned  to  move  the  probe,  the  arm  of  the 
potentiometer  rotates  a distance  proportional  to  the  probe's  motion 
along  the  axis  of  the  column.  Thus  the  voltage  of  the  arm  is  equivalent 
to  the  distance  of  the  probe  with  respect  to  the  structure.  This  voltage 
increases  continuously  as  the  receiving  structure  moves  away  from  the 
launcher. 

Both  types  of  oscillograms  were  taken  'under  different  operating 
conditions.  The  argon  discharge  was  operated  over  a range  of  neutral 

gas  pressure  and  D.  C.  discharge  current.  The  pressure  was  varied 
from  1 to  10  q.  Hg  and  the  current  between  20  and  300  mA.  In  this 
domain,  the  electron  temperature  ranged  between  2.  5 and  5 eV.  (the 
high  temperature  corresponding  to  the  lowest  pressure),  while  the 
average  electron  densities  were  of  the  order  of  10^-10^  cm"  ^ (the 
low  densities  corresponding  to  the  lowest  pressures  and  discharge  cur- 
rents). The  range  of  electron  collision  frequencies  has  been  tabulated 
in  sec.  2.  1.  1.  To  investigate  the  effect  of  column  radius  on  the  propa- 
gation of  the  waves,  column  section  of  diameters;  .6  cm,  . 952  cm 
and  1.  32  cm,  were  used. 

Finally,  the  characteristic  of  the  waves  as  a function  of  exciting 
pulse  polarity  and  amplitude  were  investigated.  A 3BG  was  also 
constructed  to  produce  negative  pulses.  For  a negative  pulse,  ions 
are  accelerated  towards  the  anode.  Since  the  amplitude  of  the  3BG 
pulse  (both  positive  and  negative)  is  constant,  attenuators  are  used  to 
obtain  the  lowest  desired  level  of  excitation  (this  corresponded  to  ap- 
proximately 100  volts).  Then,  the  peak  voltage  of  the  exiting  pulse 
(numerically  equal  to  the  longitudinal  electric  field  the  charged 
particles  between  the  plates  feel)  is  progressively  increased  (by  re- 
ducing the  attenuation)  to  its  maximum  value.  In  taking  the  oscillo- 
grams, each  time  the  exciting  pulse  is  increased,  the  received  waves 
are  attenuated  by  a corresponding  value  to  keep  the  gain  of  the  system 
fixed,  allowing  direct  comparison  of  the  profiles  of  the  received 
signals. 
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2.  2.  2 Finite  axial  magnetic  field 

To  study  the  propagation  of  large  amplitude  pulses  along  magneto 
plasma  columns,  a uniform  magnetic  field  of  strength  up  to  1,  2 K 
Gauss  was  added  to  the  set-up  of  fig.  2.  5.  The  field  was  produced  by 
seven  ESI  49  coils  arranged  in  such  a fashion  that  the  magnetic  field 
variation  over  the  test  region  was  only  _+  1%  • Moreover,  the  test 

section  in  fig.  2.  5 was  surrounded  with  a conductor,  except  at  the 
launcher,  and  receiver  locations. 

Initially,  the  same  method  (i.  e.  tube  inserted  perpendicular  to 
the  plates  through  holes  made  in  them)  was  used  to  couple  the  pulse  to 
the  magneto  plasma  column.  It  proved  to  be  very  inefficient  in  excit- 
ing propagating  bulk  waves  due  to  great  coupling  losses.  The  tube  sec- 
tion was  then  redesigned  and  copper  rings  were  inserted  such  that  the 
parallel  plates  were  in  direct  contact  with  the  plasma.  This  improved 
the  coupling  greatly  and  proved  to  be  necessary  to  launch  the  waves. 
The  modified  set  up  is  shown'  in  fig.  2.  7.  Similar  recordings  as  for 
the  case  of  no  magnetic  field  were  obtained  under  different  conditions 
of:  plasma  frequency,  magnetic  field  strength  and  background  neutral 
pressure. 

Before  presenting  the  results  of  the  experiments  (both  with  and 
without  magnetic  field),  an  analytical  formulation  and  analysis  of  the 
experimental  conditions  will  be  given  in  the  next  chapter.  In  the  light 
of  knowledge  acquired  from  such  analysis  experimental  data  will  then 
be  presented  and  analyzed. 


Fig.  2.  7 Block  diagram  of  experiment, 
finite  magnetic  field 
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III.  Analytical  formulation  of  the  experiment 

The  problem  posed  by  the  experiment  is  as  follows:  a plasma  column 
of  circular  cross-section  of  radius  b,  and  infinite  length,  either  in  free 
space  or  surrounded  by  a conductor,  is  excited  at  z = 0 by  a finite  ampli- 
tude source  of  electromagnetic  radiation.  In  general,  a finite  axial  D.  C. 
Magnetic  field  is  present  as  shown  in  Figure  3.  1.  The  spatial  and  tem- 
poral evolution  of  the  excited  bounded  modes  is  investigated.  This  analysis 
involves  the  adoption  of  an  appropriate  idealized  model  which  is  mathe- 
matically tractable. 

3.  1 Analytical  Model 

Consider  a homogeneous  plasma  column  imbedded  in  the  gas  from 
which  it  is  formed  by  partial  ionization.  The  plasma  is  considered  to  be 
a cold  electron  gas  moving  through  a stationary  neutralizing  background 
of  ions.  The  dynamical  interaction  between  electrons  and  the  background 
is  described  by  a constant  collision  frequency  v . The  interaction  is  as- 
sumed to  be  such  that  there  is  no  loss  of  electrons  due  to  ionization,  re- 
combinatibn  or  attachment.  Since  the  percentage  ionization  of  the  discharge 
is  low  (.  01%),  it  may  be  assumed  that  is  basically  the  electron  neutral 
collision  frequency.  The  fact  that  the  background  density  of  the  discharge 
is  dependent  on  radius,  undermines  the  assumption  of  a uniform  plasma 
column,  but  a compromise  is  necessary  if  the  finite  amplitude  analysis  is 
to  be  mathematically  tractable.  The  effects  introduced  by  the  inhomo- 
geneity will  be  discussed  at  the  proper  place. 

Quantitatively,  Euler's  equations  will  be  used  to  describe  the  dynam- 

(37) 

ics  of  the  electrons  : 

d fn  + V • nv  = 0 (3-1) 

mndy  + mnv  • V v + mnv  v + enE  + nv  x a H =0  (3-2) 

= electron  number  density 
= average  electron  velocity 

= collision  frequency  for  momentum  transfer  between 
electrons  and  neutrals 

m electric  field  intensity 

= background  magnetic  field  intensity 


whe  re 


n(r,  t) 
v (r,  t) 
v _ 


E(r,  t) 

H 

-.o 
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When  stated  in  the  form  of  (3-l)-(3-2),  a series  of  assumptions  must 
be  apoenaed  to  the  exact  Euler  system.  These  assumptions  may  be  found 

( O Q \ 

in  a standard  text  on  plasma  dynamics1 

Since  it  is  our  interest  to  study  the  interaction  of  the  above  plasma 
model  with  an  electromagnetic  field.  Maxwell's  equations  must  be  added 
to  (3-  l)  and  (3-2): 


x E = -u  d.  H 
- no  t ~ 


(3-3) 


v x H = e dE+J 
-v  o t ^ _e 


(3-4) 


whe  re 


= electric  current  source. 


J involves  all  electric  current  densities  present,  aonlied  or  induced.  The 
~e 

interaction  of  plasma  and  EM  field  occurs  through  the  quantities  and 
ne  E in  a self-consistent  manner.  Explicitly 


J = J + J 
_e  — app  -c 


(3-5) 


whe  re 


J 

~ app 

J 

~c 


= externally  applied  electric  current 

s electron  convection  current  in  plasma  and  is 
given  by 


J 

-c 


- env 


(3-6) 


Other  contributions  to  J (such  as  polarization,  currents,  drift  currents 


~e 


due  to  plasma  inhomogene Ltie s,  E x H drifts,  etc;  and  conduction  currents) 
will  be  neglected.  Equations  (3~l)-(3-6)  forma  determinate  system  of 
equations  through  which  any  inte raction  problem  may  be  studied,  once  J 
is  specified. 

Redefining  the  product  nv  as  a new  variable: 


u = nv 


(3-7) 


and  considering  n which  is  a function  of  space  and  time  as  composed  of 
two  parts:  a time  independent  part  or  background  and  a time  dependent 
part  or  fluctuating  term: 
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n(r,  t)  = no+  n(r,  t) 

whe  re 


(3-8) 


n 

= backgrotind  electron  density 

n 

= time  varying  electron  density 

and  n 

» n. 

o 

The 

I'JQ  } 

system  (3-  1)— (3-6)  simplifies  tok 

dfcn  + V • u = 0 

mdLu  + — y • uu=-neE-euuxH  - mu  u 
n ~ ~ ro~  ~o  ~ ~ 

o 

V x H = e d E - eu  + J 

— o t ~ - app 

V x E = - u d H 

~ o t _ 

To  this  set  of  equations,  we  also  add 

e V • E = - e(n  - n ) = - e n 
o ~ o' 

and 


(3-9) 

(3-10) 

(3-11) 

(3-12) 

(3-13) 


V • H = 0 (3-14) 

Formally  solving  equations  (3-10)  to  second  order  in  u (see  Appendix  1) 
and  using  the  results  in  eqs.  11;  the  system  is  reduced  to: 

V x H = e • E + N(E)  + J (3-1  5a) 

~ = ~ aPP 

V X E = - jUQd  H (3-1  5b) 

where  € is  a linear  operator  given  by: 

e - (€  d 1 + e H) 

= o t= 

and  N(E)  is  a nonlinear  term  in  E given  by: 

e 

N(E)  - n(  ^ V • n e iie  - — e V • E) 

m ^ c.  ~ ~ ~ ~ n - ~ 

s n ~ o 

o 

1 s unit  dyad 


with 


3- 


and 

IT  'd.H.v)  = (md  1 - eu  H x 1 - rr.v  1 ) 1 e n 

t — O U — — vj 

Decomposing  the  gradient  ooerator  as  1 - 1 ' - d z , Eq  3-15  may 

‘ z . (40) 

be  cast  into  a form  appropriate  for  g -ided  wave  analysis 
Written  in  operator  form: 

(L  + id  ru  = N u>)  + i J (3-16) 

z a pp 

where,  for  p < b (where  b is  the  radius  of  the  column) 
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( 3 - 1 6a) 


f° 
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- z x 1 
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/N(E)\ 

N ( <i  )—  i(  ] (3 -16b) 
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for  p > b,  aQ-  0 and  the  equations  reduce  to  Maxwell's  equation  in 
free  space. 

Along  with  these  defining  equations,  boundary  conditions  must  be 
specified.  These  are: 

For  a column  in  free  space 

a)  plasma  boundary  is  assumed  sharp 

b)  at  p = b,  tangential  E and  H fields  are  continuous 

c)  solutions  must  be  finite  everywhere 

d)  conditions  at  the  source  will  be  deferred  until  later 

When  the  column  is  surrounded  by  conductor,  the  above  conditions 
apply  except  that  b)  must  be  restated  as:  Tangential  E and  H fields 
are  zero  @ p = b. 

Now  we  proceed  to  obtain  a solution  to  the  above  posed  problem. 
As  a matter  of  convenience,  the  source  term  in  Eq.  3-16  will  be 
dropped  until  the  time  it  needs  to  be  considered  to  interpret  the 
experimental  results. 

Taking  Fourier  transform  in  time  of  3-16,  one  obtains: 


L (-  i u,  9.)  iM r,  oi)  + i Td  (r.io)  = N u )] 

t — Z 


(3-17) 
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where  the  transformed  operators  are  obtained  from  Equations  3-17 
by  recognizing  that  d — - - ij  in  the  transform  domain  and  that  the 
transform  of  a product  is  the  convolution  of  the  individual  transform 
(see  Appendix  1). 

<Ve  shall  seek  solutions  of  the  problem  in  the  form  of  an  expansion 

in  terms  of  the  eigenfunctions  of  the  transverse,  lossless  operator  Lg, 

where  Lq  = L | ^ _ q Assume  that  if  ( r , qj ) lies  in  the  space  spanned 
c 

by  the  transverse,  lossless  eigenvectors;  moreover,  that  the  operators 
in  Eq . ( 3-17)  act  on  this  space  . Then  if  f r , w ) may  be  represented 
as: 

if(r,w)  = ^.  y)  !fft(  r,  w) 

or 

where  the  eigenvalue  problem  for  the  is  defined: 

Lo  (-!«.%)  *a  = Ka  (u)  T 

the  eigenvectors  ijj  possess  the  orthogonality  property 

(VrV*6«p  (3-19  a) 

Using  this  property  of  the  <f  's,  the  amplitudes  in  Eq  (3-18) 
are  found  to  be  given  by: 

0^  = t4i,  Hf^)  (3-20) 

where  <f  is  the  actual  nonlinear,  lossy  field.  The  equation  for 
determining  the  co -efficients  Q may  then  be  obtained  as  follows: 

Scalar  multiply  Eq  (3-17)  by  0^  and  (3-19)  by  0 and  subtract: 

(<fa,  L0)  - (0,  + i 9Z  (r.(i,fa)  + Ka  (w)(r,fa,f ) [N(0 ), 

or 

{'i>a  <L-L0) + >■  3Z  (^.  r + /ca(<*> ) (0,  r<fa)  = [n  (0),  0J 

identifying  (L-Lq)  as  the  loss  operator  v,  using  Eq.  (3-20)  and  then 
replacing  if  (p  , u)  by  the  expansion  Eq  ~ (3-18),  we  obtain  for  the 
amplitudes  : 


(3- 18) 


(3-19) 
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i 3 G-  (z,  j ) + k (u)C'l  (z,  w ) + i7(«)ct.Jz»  J )s  i[  5.  » N(<1E)1  (3-21) 

Z Of  l<  a u U "* 

where  7 ( u ) = (>i  , 1/  1/.  ) 

Of  — Of 

Thus,  the  solution  of  the  field  problem  is  cast  into  a study  of  the 

evolution  in  z of  the  modal  amplitudes.  To  arrive  at  an  explicit  form, 

we  must  evaluate  the  two  inner  products  involved.  To  do  this,  an 

explicit  form  for  the  eigenvectors,  ^ , must  be  obtained.  Keeping 

in  mind  that  the  goal  is  in  the  analysis  of  the  experimental  results,  we 

are  interested  in  obtaining  a set  of  eigenvectors  to  represent  the  guided 

field  produced  by  an  electric  current  source  directed  along  Zq  and 

independent  of  0.  These  are  obtained  in  Appendix  2,  for  both  HQ=*  0. 

and  H =0.  The  results  are  summarized  below: 
o 

For  zero  magnetic  field,  the  discrete  spectrum  contains  a single 
eigenvector  and  the  corresponding  eigenvalue,  i.  e.  , or  = 1.  The 
properly  normalized  eigenvector,  i.  e.  , normalized  as  in  Eq.  (3-l?a), 
is  given  in  component  form: 


(Ka(u)/u<zQ)l/Z(\/2p)l/Z 


l !k  p ) 
1 A 1 


, , . .1/2,,  , ,1/2  lat\  Ibl  K1(72p) 

I ( K±  P ) 

(u  € AJ«»1/2  ( 1/2  P)^!  -u?/«2)  


p < b 


(3-22a) 


p > b 


H«e  = 


(3-22b) 


I (ft  b)  K ( k p) 

/ , ,.l/2 ,,  /,  „ ,1/2  0 1 11 

(“€o/'f«(w))  d/2,0)  K (x  b)  

0 ~ 2 X2 


i(l/«eft  /c  (w))I/2  (l/2|3)l/af  I (»q,  p)  p<b 

O Of  A 1 


( 3 - 22c) 


1 b) 

id/ we0  (w))1/5  < 1/2 p )1/2  ifrr-br  Ko(\/)p>b 

0 x2  2 
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Ilt  I^,  Kg,  K^,  are  the  modified  Bessel  functions  of  the  given  order, 

K , k and  the  eigenvalue  k(~  ) (the  a may  be  dropped  from  the 
i j a o 

eigenvalue;  however,  it  will  be  retained  on  the  eigenvector  so  as  to 
differentiate  it  from  the  total  field)  are  related  as  follows: 


2 2 2,2  2 
= *"  * (1  -u?/w  >«0 


(3-23a) 


2 2 2 

K - K - Kn 

L 1 


(1  - w “/w  2) 


*0  = u 


IAk  b)  KAk  b 

1 1 1 1 a 2 


2.  , . K K (K  b) 

K l(K  b)  i2  o a <? 

1 1 0 a i 


(3 -23b) 


(3-23c) 


For  low  Magnetic  fields  (u^  < uj^),  column  surrounded  by  conductor, 
we  make  the  assumption  that  the  waves  are  slow  and  therefore  are 
primarily  of  a TM  type.  1 The  exact  eigenvectors  are  then  approxi 
mated  by  the  quasi-static  eigenvectors: 


E = iAJ  (K  p) 
a o I , 

z M 


JAk  p) 
1 M 

K 

x 1 


A“V3 


JAK  P) 
1 M 


K 

1 1 


(3 -24a) 


(3-24b) 


(3-24c) 


Where  J , J,  are  the  Bessel  function  of  the  first  kind, 
o 1 


C3=Co 


2 2 
w -wc-> 


\ 3 1 S , K “=  - K~  — and  k b=  p 

1 M 1 €1  M n 


(3-24d) 


p is  the  nth  zero  of  J . 
n o 

Eq.  3-  19a)  by: 


The  normalization  constant  is  given  (from 
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A = 


1 


,1/2 


(3-25) 


(wKaY  e,  ) ' “ ^3bJ,  (K  b) 


In  the  limit  of  large  magnetic  fields;  i.  e.  _i  /u  > > I,  u / _j  >>1; 

c c p 

the  eigenvectors  are  obtained  from  the  exact  equations  and  are  given 
by: 


E>  = P) 

az  0 M 


( 3 -26a ) 


c.  = A 
a 

P 


k e 
a 


' P> 

* 1 i, 

M 1 


( 3 - 2 6b ) 


H«. " ’> 


9 


X1 


1 


(3-26c) 


where 


K " = e(w  (i£  - O 

ij  OO  a 


2 . 2 

€ = 1 - U)  / W 

p 


K b = o 

M 


and  the  normalization  constant  given  by: 


A’  = 


1 1 

(W/C  € )1/2  CbJ.U  b) 
a o 1 j.  l 


(3-27) 


For  completeness,  the  properly  normalized  eigenvectors  for  the  case 
of  a column  in  free  space  with  an  infinite  magnetic  field  present  are 
also  derived  in  Appendix  2,  but  not  reproduced  here. 

Using  these  eigenvectors,  the  various  terms  in  equation  3-21 
can  be  evaluated  explicitly.  Expanding  the  various  inner  product,  the 
amplitude  equation  may  be  rewritten  as: 
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id  Q (z,-j)  + k (u)<2  (z,  j)+i7(u)tf  (z,  -) 

Z Cl  U Cl 


= J jr(u,«lfu2)  ) fla(*,«2)6(u-«-  )dw  d«2  (3-28) 

-OC 

where  ^ (u  , u oj  ,?)  is  the  Kernel  of  the  integral  operator.  The 
effects  of  collisions  have  been  neglected  in  R.H.S.  y(o)  and  J?(w,  u>9) 

are  given  by  (see  Appendix  3): 

For  zero  magnetic  field,  column  in  free  space, 


2 


e<f  uj  b 

2_  f 


m u ui  ■ [p(aj»u1*w2'  p)  * w(u'“i>u2'  p)l  96(3 
1 2 o 


(3-29) 


whe  re 


P(u  , u , « 2,  p) 


K(Ml)K(u2)fC(a»)  \1/2I'1(2)  I,1(0)  I'(2) 


,3o3nlo2,,0 
^wlw2  €02  P 13  P 


«(2>  *<2>  k(1) 

M M M 


(3-29a) 


W(w,  u>lf  u2,  p)  = ^(u^) 


/C(u  ) K(u) 


\ 1 / 2 T(0)  t(l) 

__  _ \ t(2)  1 1 1 

»c(w2)«  «2we223p  plp2  / °’  fc(0)  fc(0) 

, 4 i Z O ' ij  i[ 


K(wr 


1/2 


3o3.  ,U2 

r,"-'-1“2-  '0' 


L K(w,  )k(u)oj1u9u)  e^2  p p p 


t(2) 

(1)  (2)  (0)  _1_1_  T(0) 

0 0 0 ' (2)  0 

r\ 

1 1 

(3-29b) 


f^  f(u. ) i=  0,  1,2;  and  the  prime  denotes  and 


2 

K («) 

y{u)  zii 

bl  + V 


i ,{k  pi 


,2 


n*  p) 


^ r ’ K 


pdp 


O K 


1 


For  small  axial  magnetic  fields,  column  surrounded  by  a 

conductor: 
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Ka(Ul)A(wl)  A(°j2> 


p)  p d p 


(3-30) 


this  form  of  £ is  obtained  by  using  the  assumption  that  the  maximum 
contribution  to  the  nonlinear  term  is  due  to  the  coupling  of  the  longi- 
tudinal component  of  the  total  field.  This  is  consistent  with  the  quasi  - 
static  assumption  (see  Appendix  3).  When  the  axial  field  goes  to 
infinity. 


e£ 

o 


m 


_>  to 


2 

1 


<a  (w  L ) -Al  '(oj  1 ) Al(u)  ^(gj^  ) 


p)p  dp 


(3-31) 


Note  that  equation  3-30  reduces  to  the  above  equation  in  the 

limit  w —as. 

c 

The  collision  term  is  given  by: 


- f 


A^(oj)  J“{k  p)  p d p 
0 M 


u>2i/e 
o o 

2 2 
CJ  +1/ 


(3-32) 


The  above  forms  of  X and  y apply  for  the  column  in  free  space, 
with  an  infinite  axial  magnetic  field  present.  However  the  proper 
normalization  A(w)  must  be  used  (see  Appendix  2). 

To  continue  the  analysis,  the  functional  dependence  of  on  cj 
must  be  obtained. 


3.  2 Functional  relation  between  k and  oj 

a 

a)  No  axial  magnetic  field: 

Since  there  is  only  one  mode  present,  the  functional  relationship 
may  be  obtained  from  equations  (2-23).  This  set  of  equations  is 
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transcendental  and  an  explicit  solution  for  K in  terms  of  - is 

a 

impossible.  Computer  solution  of  these  equations  is  shown  in  figure 

3.  2.  In  order  to  carry  analytical  methods  further  so  as  to  gain 

insight  into  the  phenomena,  an  approximate  analytical  function  of  w 

will  be  used  for  K^.  The  approximation  is  derived  from  the  behavior 

of  K in  eauations  2.23  for  large  and  small  os. 
a / 1 ' 

For  w — *•  0,  equation  (2-23c)  reduces  to:  ° 

99  99  99  ■’0  1/9  ^ (^P)  2 

- (^b“-  <b2)ln(«;b--  kV)1/2  - Kp^  ^ 


where  Kp  - 


« b 
P 


= — c = speed  of  light, 
c - 3 


Equation  (2 -23c)  may  be  further  reduced  to 


K 

a 


yK?> 

I ! (K?) 


1/2 


U 

UJ 

p 


(3-33) 


thus  for  oj  — 0,  k and  uj  are  linearly  related.  Moreover,  for 
a ' 

to  — “ p//T,  are  sonance  is  observed.  The  appropriate  choice  of 
approximate  equation  must  take  cognizance  of  these  facts.  Such  a 
choice  is  given  by: 


(3-34) 


where 


2 

v 


o 


2,2 
u b 

P 


2 

I (Kp)  u>  b6 

o r P 

X1(Kp)  2c2 


-1 


the  values  of  K obtained  from  (3-34)  were  compared  with  those 
obtained  by  computer  solution  of  equations  (3-23)  and  (3-34).  The 
values  of  v obtained  for  a best  fit  were: 


v 

o 


(.5163422)  - b 
P 

(.  67  165)  a b 

p 


for  14  point  fit 
for  2 1 point  fit 


(3-34a) 


u:j6 


-3 


Since  for  the  range  of  interest  equation  (3-34)  may  be  expand® 

as : 


(3-35) 


where  only  the  positive  mode  will  be  considered.  Equation  (3-35)  is 
the  desired  relationship  for  the  no  static  magnetic  field  case. 

b)  Weak  magnetic  fields;  column  surrounded  by  conductor 
From  equation  (2-24d);  solving  for  the  propagation  constant  ' 


(V 

f 

"22 
-Ui- (u“  - 

2 2* 

<J_  -<d“  ) 

P c 

1 b 

/ 

2 2 

2 2 

\ / 

f 

(u  -u“) 

(cj  -CO  ) 

L p 

c _ 

For  < u the  low  frequency  passband  of  equation  (3-36)  is 

similar  to  equation  (3-34),  except  that  resonance  now  occurs  at  u . 

c 

For  u < u < b>  and  fixing  on  modes  moving  in  direction  z > 0; 
c ? 

equation  (-3-36)  may  be  approximated  by  an  equation  similar  to  (3-35) 
as: 


= KnM  = f-  (1  + 1/2  J/J) 

s 


(3-37) 


where  v 


, u u 

_ _b p c 

“ p 2 2.1/2 

n (up+uc) 


When  the  magnetic  field  goes  to  infinity,  is  given  by: 


K 


K2=K2 ^ - 

a o r , 3/2, 

[ 1 - w / w 


again,  in  the  region  of  interest,  i.  e.  slow  waves,  the  above  becomes: 


2 

< = 


2 2,2 

/ uJ 


Of  2 / 

( 1 - U)  / Up) 
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K = K (u)  =-?-  (1+  1/  2 ) 

a n v p 


(3-38) 


Where  v =vj  / K and  ^ =o  / b.  (3 -38a 

■jo  p'  j_  j_  ■ n 

By  equations  (3-35),  (3-37)  and  (3-38),  the  coefficients  that  appear  in 
the  amplitude  equation  (3-28)  have  been  explicitly  defined.  In  the  next 
section  the  final  form  of  the  amplitude  equation  is  obtained. 
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-3*3  c inal  Form  ot  Amplitude  Equation 

Using  the  approximate  relations,  derived  in  3.2,  between  k 
and  w;  Eq.  (3-23)  may  be  written  in  its  final  form. 

A.  No  Magnetic  Field: 

Using  Eq.  (3-35),  Eq.  (3-3a),  and  the  results  of  Appendix  3 in 
Equation  (3-23),  we  finally  obtain  for  the  modal  amplitude: 

3 


id^z.u)  + (-£-  + “ 


Q n 

= C1  J Q-a(z,u2)a  (z,u  j)dx tc,  /«  2U  a u,U])a(M,)di 

o — - 1 a i (y 


O V GJ 

O D 


& jz,w)  + i -2—  a (z,«) 


vQ  a 


a 


a 


tC3  1“ 2“l  Qa(=.«2)d^  + C4  f»20i[z,»z)aa(z.»l)i>, 

(3-39) 

« 

where : d\—  5 (u  ' u 2)  du  [ dw  ^ 
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“ V,  v e 
p d • o o 


..6  3 
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v ^ 1/2 

o 
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(3- 39a) 


(3- 39b) 


(3-390 


(3-3  9d) 


and  v is  given  in  appendix  3. 

Multiplying  by  -i,  and  recognizing  that  -iu  - 3^  , the  above  equation 
transforms  to  the  time  domain  as, 
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& (z.t)  d.  & ( z , t ) - 

O'  ^ t Q 


V -J 
o p 


- d.3  Q. (2> t)  + -£■  a (2, ti  = c ar  (d.a  (z.t;)2 

Z z a v*  # I t \ t <*  / 


+ c,df(adfzfl)  +c,d?a  <*  ,t)  d 0.  iz,t)  + C Oi  (z,t)dQ,(z,t). 
2 t a t a 3 t a t a 4 a t a 

Expanding  the  R.  H.  S.  of  the  above  and  collecting  terms: 

d &(z,t)  + — d.  a (z,t)-  — -d*CL  (z,t)  + -Z£.Q.{z,t) 

z a v t a 2 t a va 

o v (J  o 

o p 

= c;  d.  a (Z,t)d2a  (z,t)+c,a  d3  a +c  a d.  a 0-40) 

ittf  t a z a t ci  -io't  ct 


where  + 


this  equation  is  the  final  form  of  the  amplitude  equation.  It  holds  for 

times  such  that  CL  (z,w)  = 0 for  w . Moreover,  the  linear  damo- 
a 

ins  term  must  be  modified  for  cases  where  Cl  (z,oj)  i 0 for  w < v 
3 a — c 

These  implications  will  be  discussed  in  section  3.4. 

B.  Infinite  Magnetic  Field. 

The  final  form  of  the  modal  amplitude  equation  for  this  case  is 
obtained  using  equations  (3-31),  (3-38)  and  the  results  of  appendix  3 
in  equation  (3.  28).  Note  that  from  Appendix  3,  it  has  been  assumed 
that  only  a single  mode  (the  lowest  order)  is  propagating  i.  e.  , the 
a 's  reduce  to  only  one.  For  the  case  of  a metal  conductor  sur- 
rounding the  plasma  column, 
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i d CL  (z,  u ) + , . -> 

z a \ v - z 

V oc  2v  “ , 

v oo  p J 


& (z,  u ) + iC  Cl  (z,  w ) 
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(k  0 ) 0 ^0 


J3(k  b) 

1 J. 


amd  Cc  given  in  Appendix  3. 
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s nice 


f j 3(k  3 


d 3 = , 72b"/2  J“  (kb) 


C',  = 1.  08  — 
4 m 


Wq 


cc 
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Jj(kb) 


multiplying  by  -i  and  transforming  into  time  domain: 

dzaaiz’t)+T~  dtQ'a(z’t)  ' T dt^a(z’t)+  Cc^z-t) 


2 V 
X 


= <='  a(z,t)dt(ZQ(z,t) 


(3-42) 


when  there  is  no  conducting  wall  present,  the  constant  above  is 
given  by: 


c r - JLjl. 


2 m 3, 

u v b 
P » 


(3-43) 


where  v^  is  given  by:  (3-38a) 

Equation  (3-28)  is  a Kortweg  -de  Vries  equation  with  a damping  term 
for  the  modal  amplitude.  It  is  valid  only  up  to  times  when  the  frequency 
spectrum  of  CL.  (z,t)  contain  frequencies  close  to  u . 


C.  Weak  Magnetic  Fieldj  column  surrounded  by  a perfect  waveguide. 

The  derivation  follows  the  same  fashion  as  that  leading  to 
equation  (3-41)  (see  Appendix  3),  with  the  exception  that  the  constant 
coefficients  are  different.  From  eqs.  (3-30),  for  the  nonlinear  term; 
equation  (3-37)  for  the  dispersion  relation  and  (3-25)  for  the  normali- 
zation constant,  we  have  from  appendix  3: 


d C (z,t)  + -—  d £*-(z,t)  - — r-  d^CL  ( 

z a v3  t or  2v  u,  2 t 


z,  t) 
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= C''CL  (z,  t)  d.  Cl.  (z,  t) 

4 O'  t at 


(3-44) 
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where  Cc  a collision  coefficient  and  is  approximately  given 
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3.  4 Analysis  of  the  Amplitude  Equation 

Comparing  the  magnitude  of  the  terms  on  the  right  hand  side  of 
equation  (3-39).  we  find  that  the  last  term  [the  term  whose  coefficient 
*-s  t*le  dominant  one.  From  equation  6-37of  Appendix  3,  ob- 

serve that  this  term  arises  from  the  nonlinear  coupling  of  the  longitu- 
dinal electric  field  of  waves  of  different  frequencies.  Since,  from 
equations  (2-2  1),  the  magnitude  of  the  longitudinal  field  is  greater  than 
the  transverse  field,  it  is  expected  that  the  strongest  nonlinearity  would 
result  from  interactions  involving  only  the  longitudinal  components  of 
the  field.  Neglecting  the  other  two  terms  in  the  R.  H.  S.  of  (3-39),  the 
resulting  equation  has  the  same  form  as  equation  (3-42),  i.  e.  , Eq. 
(3-39)  reduces  to  the  K-de  V equation. 

Thus  the  evolution  in  z of  the  modal  amplitudes  is  described  by 
an  equation  of  the  K-de  V type.  This  results  in  the  recognition  that 
propagation  of  large  amplitude  slow  waves  and  surface  waves  along  a 
plasma  column  belong  to  a very  general  class  of  nonlinear  wave  phe- 
nomena, i.  e. , nonlinear  dispersive  waves.  It  is  important  to  note 
that  the  sign  of  the  nonlinear  term  in  (3-42)  is  opposite  to  the  corres- 
ponding one  in  (3-39).  This  difference  results  in  qualitatively  different 
behavior  of  the  solution  which  will  be  discussed  in  the  next  section. 

Noting  from  Equation  (3-18)  that,  since  for  each  case  we  are 
only  considering  a single  mode,  the  field  solution  is: 

w 

* (P.  t;z)  = f ajz,  W)  </,  (p,  w)^4 

J Ce  G £7 T 

-oo 

where  z,  w)  and  < //  (p,w)  are  given  by  Eqs.  (3-39)  and  (3-22)  res- 
pectively, for  surface  waves  and  by  Eqs.  (3-24)  and  (3-42)  respective- 
ly for  body  waves.  In  particular,  for  the  longitudinal  electric  field 
of  the  wave  s : 

E (o,  u»;z)  = 2?  (z,  w)  E (p,  w) 

^ Cc  a 

Putting  in  the  value  for  E^(p,  u>)  for  either  case,  we  find  in  the  regime 
Kkb  < 1 


E (p,  u;z)~  const  • iu  &(z,  u) 

z Ct 


(3-45a) 


-39- 


or  from  inversion: 

E (p,  t ;z)  „ - const  3t  (z>  t). 


(3 -45b) 


This  result  is  very  convenient  since  it  allows  the  solution  of  Eqs. 
(3-42)  or  (3-39)  to  be  immediately  associated  with  the  longitudinal  elec 
trie  field  of  the  wave.  Thus  the  Modal  Amplitude  acts  as  a potential 
whose  time  derivative  is  proportional  to  the  longitudinal  E field. 
Rewriting  equations  (3-42)  and  (3-39)  as 

az  GLa  (z,  t)  + -i-  3t  (z,  t)  - cd  3 t &a  (z*  ^ ~ Cc  ^ct  ^Z’ ^ 


-c"  t)  at^,t)=o 


(3-46) 


whe  re 
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With  the  observation  made  in  (3-45),  Equation  (3-46)  may  be 

used  as  the  starting  point  for  the  study  of  pulse  propagation,  of  limited 

band-width  and  finite  amplitude,  along  a plasma  column  with  infinite 

axial  magnetic  field,  or  with  a weak  magnetic  field,  respectively, 

34 

Dcezey  et  al  , using  quasi-static  approximation  derived  a K-de-V 
Equation  for  the  evolution  in  time  of  the  potential  of  an  electron  acous- 
tic wave  propagating  in  a cylindrical  plasma  column  in  an  infinite  mag- 
netic field  and  surrounded  by  a perfect  conductor. 

The  K-de-V  equation  can  be  obtained  for  other  kinds  of  plasma 
waves,  for  example  ion  waves  in  a cold  plasma  for  the  case  w 1^} 

hydromagnetic  waves^^,  and  for  other  wave  phenomena^2-,  For  a 
comprehensive  survey  of  K-de-V  literature  see  [44]. 

As  previously  pointed  out  we  could  not  operate  the  experiment 

in  the  regime  u » u Results  were  obtained  only  for  u =0  and 
ce  pe  ’ ce 

< w . Even  though  the  dispersion  relation  for  symmetric  body 
waves  in  the  regime  fce  < f can  be  approximated  by  an  equation  of 
the  form  (3-37),  the  structure  of  the  properly  normalized  and  exact 
eigenvectors  is  rather  complex  and  is  not  amenable  to  analytical  mani- 
pulation. The  coupling  coefficients  in  the  equations  for  the  modal 
amplitudes  (more  than  one  mode  has  to  be  conside  red^  ^ ) would  not  be 
simple  functions  of  frequency,  and  numerical  solution  would  have  been 
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inevitable.  By  properly  approximating  the  exact  eigenvectors,  the 
coupling  coefficients  in  the  modal  amplitude  equations  could  be 
simplified,  and  the  analytical  result  in  the  form  of  Equation  (3-4oi 
was  then,  possible. 

Since  no  attempt  is  made  here  to  solve  the  K-de-V  equation 
either  analytically  or  numerically  the  equation  for  the  modal  ampli- 
tude will  be  subsequently  normalized  30  that  the  results  obtained  by 
others^-"^  may  be  used  to  enable  us  to  rephrase  the  results  obtained 
in  Chapter  IIL 

Introducing  new  variables,  and  restricting  ourselves  to  experi- 
mental conditions: 


^ (2,  t)  = 

or 


4-  uipb  I I dQ  ( 2,  t)  No  magnetic  field 


Jcb!  C4  I Cla  (z> 


Weak  magnetic  field 


(3-47a) 


z'  = 


t1  = 
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*p(t  - z/vq) 
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No  magnetic  field 
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No  magnetic  field 


Weak  magnetic  field 


We  obtain  from  (3-46)  and  (3-47): 

az,  ^(*'(t')+bc^(*',f)  dt,aa (z'.f) 

- or'8  J fijz1,  t')  = 0 


*(3-47b) 


(3-47c) 


(3-48) 


No  magnetic  field 


whe  re  or'  = 


Weak  magnetic  field 
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Equation  (3-48)  is  in  the  "standard  form"  (neglecting  the  collision  term) 
discussed  by  Berezin  and  Karpman1  ' J.  They  investigated  the  evolu- 
tion, tor  different  values  of  a , of  a disturbance  at  z - 0,  radiated  for  a 

/\ 

bounded  time  interval  A , i.  e. , given  <2.  (0,  t' ) ~ & (V  / A).  How  does 

Ct 

32^(2'.  t')  depend  on  A and  o',  if  Cl  (z'f  t' ) evolves  as  in  (3-43). 

The  boundary  condition  for  (3-48)  may  be  obtained  from  the  ex- 
periment as  follows:  from  figure  2-7,  as  the  BBG  pulse  propagates 
in  the  parallel  plate  structure,  part  of  its  energy  is  coupled  to  a radia- 
tion field,  part  to  a plasma  guided  field  and  the  rest  is  dissipated  by 
the  termination.  The  amount  of  energy  that  goes  into  each  type  depends 
on  how  well  the  parallel  plate  couples  to  that  mode.  Even  though  we 
measured  the  coupling  of  energy  by  the  parallel  plate  structure  to  the 
guided  field  as  a function  of  frequency,  an  absolute  measurement  on 
the  percentage  of  energy  of  the  initial  pulse  that  goes  into  the  guided 
field  could  not  be  done.  The  unnormalized  modal  amplitude  at  z = 0, 
i.  e. , <fZ^(0,  t)  is  related  to  the  longitudinal  electric  field  component  of 
guided  field,  excited  by  the  pulse,  by  equation  (3-45a),  and  therefore 
from  equation  (3-45b): 

a)  No  magnetic  field: 


3t.  d.J0,f)  = 


1_ 

2 


(b,  0,  t1  ) 


(3-49) 


Note  that  the  factor  in  parenthesis  has  the  dimensions  of  (coul-sec^)/ 
(kg-cm)  or  of  (electric  field)"^  as  it  should.  Since 


v = * b(.  67165): 

o p 


3t,^?(0,  t')  = - 1.  9503  x 10"  E (b,  0,  f 

t Ct  Cm-  Z 
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E in 
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Coul-  sec' 


(3-50) 
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E (b,  0,  t' ) is  the  longitudinal  field,  evaluated  at  the  boundary  of  the 
column  of  the  "already  set  up"  guided  field  at  z - 0.  It  is  not  the  ac- 
tual B3G  pulse.  However,  since  the  field  is  applied  for  a finite  amount 
of  time,  we  expect  that  the  guided  field  at  z = 0 also  possess  this  char- 
acteristic in  time,  but  as  mentioned  before,  its  absolute  strength  is 
unknown. 

b)  Weak  magnetic  field 

From  (3-24),  E (p,  0,  w)  ~ E (0,  oo)  J (k  p) 

Z 7.  O j_ 

^ ,1  + J. 

(0,  t')  = + (3.  664  x 1 0")  -B--5  --  EJ0,  t')  (3-51) 

z ct  c z.  , z 

cj  uj  b 
p c 

Looking  into  the  physical  consequences  of  the  above  theory,  let  us 
first  consider  the  linearized  equation  obtained  from  (3-4 3)  when  the 

/N  /V 

nonlinear  term,  CLa  , is  neglected.  The  resulting  equation  is 

the  Airy  equation  with  a damping  term 

dz,  t')  +bc'cz£(z\  f)  - aa(  z',  t')  = 0 


whose  solution  is: 

a (z',t')=  (S  («)  ei(/d“)z'  ■“t')  du 
a M a 17 

with 


(3-52) 


k(w)  = - ib  C 


and  the  normalized  amplitude  A (u)  is  given  in  terms  of  the  boundary 
condition  (3-50)  or  (3-51): 

30 


Vw)  = lio  I (const)  Ez(0it)  e 


iut 


dt 


-00 


and  from  (3-45),  we  get  for  the  actual  guided  field: 


E (z\  f)  = f (const)  E (0,  3j)  e1^^2'  "ut' > p- 
2 .io  2 2:r 


or 
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E (z'.f  ) = (const)  e'DCcz'  J G(z',  t1  -t" ) E (0,  t')dt" 

-00 


(3-53) 


with 


G(z\  f)  = 


Ai  - 


\ 


(3^z»)1/3  ^ (, 2a'z')l/ 3 


wnere 


Ai(x)  = ± f ei(r3/3  +XT: 
-00 


dT 


Since  the  initial  width  of  the  applied  signal  is  sufficiently  small, 
then  for  large  z'  , G(z',  t1 ) is  a slow  varying  function,  and  the  solu- 
tion (3-53)  of  the  linearized  equation  can  be  written  in  the  following 
form'42': 


-bC  2' 
c 


Ez(z' , t1 ) ~ (const) 


Ai 


Oa’z'  )1/3  V (3^z*  )1/3J 


for  large  z'(3-54) 


E^  does  not  depend  on  the  character  of  the  initial  radiation  as  long  as 
its  width  is  small.  However,  for  small  z'  , the  solution  depends 
strongly  on  the  boundary  condition.  The  solution  (3-54)  represents  a 
wavepacket  moving  away  from  the  launcher  with  the  first  peak  moving 
with  velocity  vq  for  Bq  = 0,  and  v^  for  0.  The  packet  spreads 

and  its  amplitude  is  reduced  as  it  moves.  Figure  3-3  shows  a plot  of 
the  velocity  of  the  first  peak  as  a function  of  plasma  parameters,  i.  e.  , 
ujp  and  u)c<  The  solution  (3-54)  is  discussed  further  in  Chapter  III  in 
relation  to  the  experimental  results. 

Next,  for  the  full  nonlinear  equation,  the  nonlinear  narameter 

r 42 1 

<t  ■ J provides  a measure  for  the  nonline  ar  effects.  To  compute  cr  for 

typical  experimental  conditions,  assume  the  boundary  conditions  for 

Eq.  (3-48)  to  be  given  at  z = 0 + Az,  i.  e. , a few  centimeters  to  the 

right  of  the  launching  plates.  Placing  the  receiving  plates  close  to  the 

launcher,  we  observe  a wave  packet  of  four  oscillations  extended  over 

an  interval  of  time  r . Thus,  take  as  the  boundary  condition  for 
— \ 

/7  (z‘,  t1 ) a function  of  the  tvpe: 

Ct 


CLQ  (Az1,  t’)  = + cp{ t1  /tw.) 


(3-65) 


1 st  peak  velm-  ity  (10 
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A1  and  31  are  unknown  constant  amplitudes  in  volts/cm  which  depend 
on  how  well  the  BBG  pulse  couples  to  the  mode  field.  From  [42], 

<T  = 1 

V? 

Putting  in  for  & and  or1: 

For  B = 0; 

_o 

it  = 4.  455  x107t  -/  A1’  (3-5” 


For  ?o\°: 


' / Z , 2 

cr  = 7.  4467  x 10  t ( — ^ 

V*'“p 


vrBr 


Taking  widths  of  the  order  of  5 nonsec.,  w ~ 10.  49  x 10~  rad/sec.  , 

9 . P 

~ 5.  18  x 10  rad/ sec  and  a ” L.  4 (see  figure  3-4  for  the  behavior 
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From  numerical  solutions  of  the  K-de-V  equation,  it  is  known"  ' that 
a critical  <r  = <r  = i/T?  exists  in  the  sense  that  qualitatively  different 

C UZl 

solutions  are  obtained  for  <r  » cr  and  <r  <<  t l j.  For  j*  > cr  , 
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nonlinear  stationary  solutions  of  the  "solitary  type"  are  permissible, 
whereas  for  t < .r  "nonsolitary"  solutions  corresoonding  to  raoidly- 
oscillating  wave  packets  are  obtained. 

Experimentally,  the  maximum  field  amplitude  available  at  the 
plates  (due  to  the  4.  4 dB  attenuation  of  the  60  nansec.  delay  line)  is 
1.  92  kilovolts/cm.  Even  though  absolute  field  amplitude  of  the  excited 
wave  is  unknown,  i.  e. , A'  or  B'  , the  first  nonlinear  effects  which 
begin  to  be  observed  is  in  the  region  cr  < r Therefore  in  the  first 
assessment  of  nonlinear  effects,  we  proceed  as  in  [45]  by  expanding 
the  solution  of  the  K-de-V  and  the  wave  number  ast 


^(z,  t)  = e CL^Q  ) + €~  <^(9  )+. . . 

k (u>)  = (u)  + e /c,  (w)  + e “ k,(w) +. . . 

at  o i c. 


(3-60) 


where 

6 = K z - wt  and  e = € (<f?  ). 
ur  a 

and  /c-  (uj ) is  chosen  to  be  zero  to  avoid  a secular  exnression  of  Cl  . 
l a 

In  a straightforward  fashion,  we  obtain  for  the  coefficients  of 
the  expansion: 

- cos  9 Cl-,- j cos  Id 

"12  Cd  uj  bu 
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K (u)  - — a+Cda 

C V. 
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(3-o  1) 
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24  Cdu  ^ b^  w 
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The  critical  difference  between  these  results  and  the  linear  results  in 
equation  (3-52)  is  the  dependence  of  the  dispersion  equation  on  the  am- 
plitude. The  specific  consequence  of  this  observation  is  that  the  char- 

f 45] 

acteristics  of  the  modulation  equations1  J become  doubly  degenerate, 
i.  e. , it  predicts  the  eventual  splitting  of  a wave  packet  into  two  separate 
ones  propagating  with  characte ristic  velocities  given  by: 
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ror  the  slow  wavepacket 


«o(w)  -T7 


*oM  - TT ~S 

p 


(3-62) 


for  the  fast  wavepacket. 


where  CL/  is  the  normalized  amplitude  of  the  nonuniform  wavetrain,  and 
o 

k'  (w)  = K . Using  these  velocities,  we  can  get  an  upper  estimate  as  to  th 
o'  dtu  o * n 

distance  away  from  the  source  at  which  this  separation  can  be  observed:1 


C+-C- 


or  using  (3-62): 

cjj  b t 
c 

2 = 

-2  a 


For  weakly  magnetic  field,  from  equations  (3-47)  and  (3-56a)!  Note 
that  for  surface  waves  we  obtain  a negative  value  for  z.  The  conse- 
quence of  this  result  is  discussed  in  the  next  chapter.  ) 


3 K“ 

u>  b 

c 


3.66  xlO15  u } + w2 
P c 


For  the  typical  experimental  parameter  previously  considered! 
2.  06  x IQ3 


(3-63) 


thus  for  a wavepacket  of  actual  field  strength  of  50  volts /cm,  separa- 
tion of  approximately  3 nansec.  will  result  at  61  cm  away  from  the 
bounda  ry. 

These  results  are  used  in  the  analysis  of  the  experimental  ob- 
servations which  are  discussed  next. 


Ul  (l 
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IV.  Results 

In.  light  of  the  analytical  results  derived  in  chapter  III,  the  data  ob- 
tained using  the  apparatus  described  in  chapter  II  will  be  subsequently 
presented  and  analyzed.  The  results  for  the  no  anial  magnetic  field  case 
are  presented  in  section  4.  1;  and  in  section  4.  2,  the  results  for  the  weak 
axial  magnetic  field  case  are  presented.  As  previously  mentioned,  there 
are  two  fundemental  differences  between  these  two  cases: 

1.  Even  though  a Kortweg-de  Vries  type  equation  describes  the 
evolution  in  space-time  of  a given  disturbance,  different  qualitative  re- 
sults will  be  obtained  in  both  cases  for  the  same  excitation.  This  is  due 
to  the  fact  that  the  signs  of  the  nonlinear  terms  (or  of  the  normalized 
initial  distribution  in  Eq.  3-46)  are  opposite.  In  interpreting  the  results 
of  this  chapter,  it  is  important  to  note  that  the  Electric  field  at  the  launch- 
ing plate  points  in  the  - .z  direction  for  what  has  been  referred  to  as  the 
positive  pulse  in  the  preciding  chapters. 

2.  The  coupling  structures  used  for  either  case  were  different. 

The  set  up  shown  in  Fig.  2-5  was  used  for  the  no  magnetic  field  case, 
where  as  for  the  weak  magnetic  field,  the  set  up  is  shown  in  Fig.  2-7. 

In  the  latter  case,  the  addition  of  the  brass  rings  allowed  for  the  pulsed 
field  to  be  directly  coupled  to  the  plasma,  improving  the  launching  ef- 
ficiency for  the  body  waves.  For  Bo  = O,  surface  waves  were  strongly 
excited  without  the  need  of  the  brass  rings.  For  these  waves,  the  fields 
are  concentrated  at  the  plasma  boundary. 

Each  of  the  above  sections,  i.e,  , 4.  1 and  4.  2,  are  divided  into  two 

subsections.  The  low  amplitude  or  linear  results  are  presented  in  the 
first  subsection;  while  in  the  second  subsection,  we  present  the  finite 
amplitude  or  nonlinear  results. 

4.  1 Zero  Axial  Magnetic  Field 
4.  1.  1 Linear  Regime 

For  experiments  in  this  regime,  the  peak  pulse  electric  field  at  the 
launching  plates  was  kept  below  192  Volts/ cm.  At  these  field  levels,  the 
integral  solution,  equation  (3-53),  obtained  in  section  3.4  is  valid.  From 
this  equation,  the  mechanism  that  affects  the  propagation  of  the  pulse  is 
dispersion. 

The  vital  experimentally  controlled  parameters  which  appear  in  the 
solution  (3-53)  are:  the  background  neutral  pressure  (i.e.  collision  fre- 
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quency),  initial  pulse  polarity,  plasma  frequency  and  plasma  column 
radius.  As  a matter  of  organization,  the  effect  of  each  of  the  above 
parameters  on  the  propagation  of  the  pulse  will  be  sequentially  pre- 
sented and  quantitatively  compared  to  the  solution  (3-53). 

Pressure  effects  had  no  particular  quantitative  interest  to  us  ex- 
cept from  the  point  of  view  of  being  able  to  launch  the  waves.  Due  to 
the  increase  with  pressure  of  electron-neutral  collision,  it  was  difficult 
to  launch  the  waves  as  illustrated  in  the  sequence  of  Fig.  4.  1.  High 
attenuation,  column  striations,  and  loss  of  statistical  sampling  coherence 
are  the  factors  responsible  for  this  behavior.  For  the  rest  of  the  ex- 
periments, the  pressure  was  kept  at  a point  where  collision  related  ef- 
fects had  the  least  influence  on  the  wavepacket  characteristics. 

The  effect,  in  this  regime,  of  pulse  polarity  is  trivial:  it  merely 
changes  the  sign  of  the  solution  (3-53).  In  Fig.  (4.2a),  we  display  os- 
cillogram in  time  obtained  at  z = 56  cm  for  the  "positive"  pulse  and  in  Fig. 
(4.  2b)  for  the  "negative"  pulse.  Figure  4.2b  is  just  the  negative  of  Fig. 

4.  2a.  The  polarity  of  the  first  oscillation  is  as  predicted  from  the  solution 
(3-54).  Quantitative  differences  rn  the  freguency  of  oscillation  and  am- 
plitude are  due  to  slight  differences  in  plasma  frequency  and  initial  pulse 
characteristics.  In  Fig.  a,  the  "positive"  pulse  generator  produces,  at 
the  launching  plates,  an  Electric  field  pulse  in  the  negative  z direction, 
its  strength  is  approximately  192  volts/cm  and  its  duration  is  approxi- 
mately 1 nanosecond.  The  "negative"  generator,  in  the  other  hand,  pro- 
duces a field  in  the  positive  z direction  and  its  duration  is  approximately 
. 5 nanosecond.  The  seemingly  reverse  naming  comes  from  the  fact  that 
the  labeling  is  appropriate  for  the  voltage  of  the  front  plate  with  respect 
to  the  ground  plate.  The  evolution  in  space  of  the  initial  disturbance  is 
shown  in  Fig.  4.  3.  The  a)sequence  corresponds  to  the  positive  pulse 
whereas  the  b)sequence  correspond  to  the  negative  pulse.  This  sequence 
is  explained,  in  terms  of  equation  (3-5  3)  as  follows:  at  z = o,  the 
launcher  position,  a wave  is  radiated  for  a bounded  time  interval.  As  the 
pul3e  travels  to  the  right,  i.e.,  z > 0,  it  rapidly  breaks  into  oscillations 
due  to  strong  dispersive  effects.  Since  the  dispersion  is  negative,  the 
long  wavelength  are  observed  to  propagate  at  a higher  speed. 
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Figures  4.4  a,b,  and  c show  oscillograms  in  time  for  different 
positions  of  the  receiver.  The  part  of  the  oscillogram  marked  with  an 
A corresponds  to  free  space  electromagnetic  radiation  propagating  at 
the  speed  of  light.  This  radiation  is  observed  even  if  the  plasma  were 
not  present  (see  Fig.  4.2  ).  From  these  curves,  we  again  observe  the 
negative  dispessive  properties  of  the  system.  This  property  is  further 
displayed  by  fourier  analyzing  the  wave  forms  of  Fig.  4.4.  From  the 
relation 

It  Aw  ~ 1 

where  At  = spread  of  the  wavepacket  in  time 

Aw  a spread  in  frequency 

Thus,  as  the  packet  propagates  and  spreads  in  time  due  to  dispersion, 
its  fourier  transform  becomes  narrower.  This  is  displayed  in  Fig. 

(4.  5). 

As  the  averaged  plasma  frequency  is  changed,  the  packet  charac- 
teristics, are  also  observed  to  change.  The  speed  of  propagation,  fre- 
quency of  oscillation,  width  of  the  packet  and  the  normalized  amplitude 
increase  as  a function  of  f^.  Time  oscillograms  for  various  averaged 
plasma  frequencies  are  shown  in  Figure  (4.6).  Figure  (4.7)  shows  a 
plot  of  the  velocity  of  the  first  peak  of  the  wavepacket  as  a function  of 
plasma  frequency.  Two  different  methods  were  used  to  calculate  this 
velocity.  The  first  uses  the  fact  that  the  first  observed  disturbance 
corresponds  to  the  free  space  electromagnetic  radiation  travelling  at  the 
speed  of  light.  By  comparing  the  time  of  arrival  at  the  point  of  observa- 
tion, of  the  first  peak  of  the  wavepacket  with  respect  to  the  EM  signal, 
the  speed  of  propagation  is  obtained.  The  second  method  uses  the  time 
of  flight  of  the  first  peak  between  two  points.  Also  plotted  in  Figure 
(4.7)  is  equation  (3- 34a).  As  was  pointed  out  in  section  3.4,  the  speed 
of  propagation  of  the  first  peak  is  given  by  v . The  experimental  values 
differ  from  the  theoretical  by  20%.  This  difference  lies  in  the  manner 
that  vq  was  obtained,  i.e.  by  numerically  fitting  the  solution  of  the  dis- 
persion relation  to  equation  (3-34a).  In  Equation  (3-34a),  14  points  were 
used  to  determine  the  low  frequency  or  linear  slope  of  equation  (3- 34a), 
i.e.,  vq.  If  fewer  points  are  used,  a large  vq  is  obtained  (see  equation 
(3- 34b)  for  example)  and  a better  approximation  to  the  propagation  speed. 
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For  receiver  positions  far  away  from  the  launcher,  i.e.,  at  large 

values  of  z,  the  far  field  solution,  equation  (3-54)  is  applicable.  This 

imolies  that  the  observed  oscillograms  should  benave  as  Airy  functions 

1/3 

of  argument  u ( t - z/  v ) / ( 3 a-1  z/  b)  . Using  the  asymptotic  ap- 
spo 

proximation  for  the  Airy  function,  the  solution  (3-54)  takes  the  form: 


E 

z 


1 sin  2/3  t'3/2  4 n/4jt>z/vc 


where  t'  = « (t-z/v  ) / ( 3or‘z  ' ) ' z'-z/b 
P 0 

the  zeros  of  the  above  equation  occur  at: 

4'  = (n  - 1/4)  2/ 3 n^3(3»2')^3  (4-1) 

These  values  can  be  directly  compared  to  the  properly  normalized  zeros 
of  the  observed  oscillograms.  In  Table  (4-1),  we  tabulate  the  zeros  ob- 
tained from  equation  (4-1)  and  from  two  different  oscillograms.  These 
are  also  plotted  in  Figure  (4-8).  If  instead  of  using  the  average  plasma 
frequency  as  measured  using  the  Cavity  method,  a plasma  frequency  is 
chosen  as  to  fit  the  nth  experimental  crossing  to  the  nth  theoretical 
crossing;  we  observe  that  all  other  (n-1)  crossings  fall  to  within  6%  of 
the  theoretically  predicted  values.  Larger  'deviations  are  obtained  for 
the  a = 0 and  n = 1 crossings  where  he  asymptotic  expansion  does  not 
orovide  an  accurate  representation.  The  value  of  « obtained  in  this 
fashion  is  20%  off  from  the  value  measured  by  the  Cavity  method.  The 
discrepancy  lies  in  the  fact  that,  the  Cavity  method  measures  the  aver- 
age plasma  frequency  of  a non  uniform  plasma  column,  whereas  the  +>  ^ 
obtained  by  matching  the  nth  zero  is  obtained  from  a theory  hat  assumes 
an  uniform  column.  The  plasma  frequency  thus  derived  corresponds  to 
a value  near  he  edge  of  the  non  uniform  column.  Thus,  we  have  de- 
veloped a method  for  measuring  he  plasma  frequency,  near  he  edge  of 
he  column,  from  he  zero  crossings  of  he  observed  time  oscillograms. 
It  has  a higher  resolution  than  he  Cavity  method  as  seen  from  Figs. 
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Table  4.  1 


Parameters: 


Receiver  distance  (z) 
Column  Radius  (b) 

Normalized  distance  (z') 

/ 

a 

v^  (experimental) 


54.  375  cm 
. 66  cm 


82.  386 
1. 48887 

Q . 

3.395x10  cm/  sec. 


f (from  cavity)  = .9 GHz 

P 


Zero  crossing  (n) 

Experimental  t' 
(u)_  measured 
using  cavity) 

Experimental  t' 

(u>_  chosen  so  as  to 
fft  nth  zero  crossing) 

j 

(a) 

(b) 

0 

9. 14898 

12.  16 

9.  623 

9.  8 

1 

19.02 

23.  16 

22.  23 

22.7 

9 

33. 43 

43. 34 

34.  33 

34.  99 

3 

45.  19 

57.11 

* 45.  19 

46.  0 

4 

55.  56 

69.  33 

54.  35 

55.9 

5 

65. 020 

30. 

63.  30 

64. 528 

6 

74.  01 

90.  65 

71.71 

73.  10 

7 

82. 3131 

102. 07 

80.  75 

* 82. 3131 

8 

90. 2908 

111. 23 

88. 

89.7 

9 

97.399 

121.1 

95.  82 

97.68 

10 

105. 1395 

130. 293 

103. 08 

105.06 

11 

112. 23 

139. 4490 

110. 32 

112. 45 

119. 158 

1 48 . 55 

1 17. 52 

119.  79 

125.821 

156. 187 

123. 56 

125.  9 

nth  zero  matched 


(a)  = .712GHz 

f corresponding  to  column  1 


(b)  = . 725  GHz 


Normalized  'Lime  position  of  zeros  of 
Electric  Field 


Figu  re 
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(4.3).  The  average  plasma  frequency  as  measured  by  the  Cavity -was 
.9  GHz  for  both  oscillograms,  whereas  due  to  small  background  pres- 
sure change,  the  number  density  had  drifted  slightly  causing  a change 
in  the  propagation  characte  ristics  of  the  wavepacket.  The  drift  cor- 
responded to  a 5%  change  which  was  not  detected  by  the  Cavity. 

Qualitatively,  the  behavior  of  the  envelope  of  the  observed  oscillo- 
grams is  as  predicted  by  equation  (3-54).  However,  quantitatively,  the 
rates  of  decay  of  the  envelopes  do  not  agree.  In  Figure  (4-9),  we  dis- 
play the  envelope  amplitude  as  a function  of  time.  Two  rates  of  decay 
are  observed.  For  rates  just  to  the  right  of  the  maximum  amplitude, 
the  decay  goes  as  t ^ ^ , whereas  for  large  times  the  rate  goes  as  t * . 

The  asymptotic  result  of  equation  (3-54)  predicts  an  uniform  decay  rate 
_ i / q 

of  t ' To  account  for  the  difference,  two  fundamental  assumptions 
made  in  the  derivation  of  the  theory  must  be  modified  if  a quantitative 
description  of  the  amplitude  is  desired.  First,  in  the  desivation  of  the 
collision  term,  the  dependence  of  (c'tf>  v if'a)  on  frequency  was  neglected 
in  order  to  get  at  analytical  results.  Secondly,  sheath  effects  which  in- 
troduce strong  landau  dampling  of  the  slow  waves  at  high  frequency 
must  be  included  in  the  theory.  The  high  attenuation  ( — t *)  observed 
for  the  high  frequencies  is  an  indications  that  this  mechanism  may  be 
present. 

The  last  parameter  considered  was  tube  radius.  Figures  4-1  Oa 
and  4- 10b  display  recordings  in  time  and  space  obtained  using  a tube  of 
radius  . 32  5 cm.  The  behavior  of  the  packet  is  as  predicted  by  the  theory 
of  chaoter  III  with  the  proper  value  of  radius  b.  A large  discrepancy 
is  noted  (see  Fig.  4-11)  between  the  zero  crossings  of  the  observed 
oscillations  normalized  to  the  Cavity  measured  plasma  frequency  and 
those  computed  from  equation  4-1.  Also  shown  in  Fig.  4-11  are  the 
crossings  computed  by  normalizing  the  experimental  nth  crossing  to  the 
corresponding  theoretical  value.  The  large  difference  for  this  tube  be- 
tween the  Cavity  measured'*)  and  the  nth  crossing  computed*)  is  con- 

P P 

sistent  with  the  previous  explanation  given  for  this  difference.  It  lies  in 
the  fact  that  the  small  column  has  a stronger  transverse  inhomogeneity 
in  the  number  density  than  the  larger  tube. 
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a)  f = .79  GHz  z * 56  cm 

? 


b)  f * .79  GHz  z - 27  cm. 
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Fig.  4.10  Oscillograms  for  differen  plasma  frequencies.  Tube  radius 
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4.1.2  Nonlinear  Regimes 

As  the  amplitude  of  the  initial  pulse  in 
start  to  play  a role.  Experiments  in  this  regime  were  carried  out  using 
peak  field  strengths  greater  than  600  volts/ cm,  at  the  launching  plates. 

It  is  important  to  note  that  for  all  experiments  (both  zero  and  finite  mag- 
netic fields),  the  signals  observed  in  the  oscilloscope  are  normalized  so 
that  if  the  response  were  linear,  the  observed  oscillogram  as  a function 
of  input  pulse  strength  will  coincide.  This  is  accomplished  by  keeping 
the  total  attenuation,  in  the  path  the  signal  travel,  constant.  Thus  to  go 
from  the  linear  to  the  nonlinear  regime,  the  attenuation  is  transferred 
from  the  input  of  the  launching  plates  to  the  input  of  the  oscilloscope.  In 
this  fashion,  effects  that  are  clearly  amplitude  dependent  are  identified. 

In  Figure  4-12,  the  profile  in  time  of  the  wavepacket  is  shown  as 
a function  of  input  pulse  amplitude.  The  positive  pulse  was  used  as  the 
source.  Qualitatively,  these  diagrams  are  similar,  although  quantita  - 
tively  they  differ  in  the  law  governing  the  decrease  of  the  amplitude  in 
time  and  in  3pace.  Moreover,  the  characteristics  of  the  first  two 
oscillations  have  changed.’  Note  that  the  maximum  normalized  ampli- 
tude of  the  nonlinear  response  is  less  than  for  the  linear  case,  implying 
that  the  surface  wavepacket  has  reached  saturation.  In  Figures  4-13  and 
4-14,  the  behavior  of  the  amplitude  as  a function  of  space  and  time  for 
the  linear  and  nonlinear  cases  are  displayed.  The  slower  decay  rate  ob- 
served for  the  non  linear  case  is  explained  from  the  dynamics  of  the  pro- 
pagation as  described  by  the  model  developed  in  section  3.  3.  For  posi- 
tive oulse,  E (0,  t) 1 E (0,  t)  i,  where  E is  the  pulse  field  strength. 

z p P 

Thus  for  any  value  of  <r , the  solution  of  the  K de  V equation  is  always 
oscillatory  and  the  character  of  the  solution  is  very  similar  to  the  linear 

T 4?  1 

result.  L “ J A physical  understanding  of  the  processes  taking  part  in 
the  propagation  of  the  high  amplitude  pulse  may  be  gained  from  the 
fourier  transform  of  the  observed  oscillograms  as  a function  of  plasma 
frequency  and  receiver  location.  In  Figure  4-15,  we  display  side  by 
side  the  nonlinear  response,  and  the  corresponding  power  spectrum 
as  a function  of  plasma  frequency.  In  Figure  4-16,  the  plasma  frequency 
is  kept  constant  and  the  recordings  are  shown  for  3 different  positions. 

These  diagrams  show  the  channeling  of  energy  from  the  lower  fre- 
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quencies  to  the  higher  frequencies  which  in  turn  are  strongly  attenuated  as 
explained  in  section  3-3.  This  accounts  for  the  saturation  phenomena  ob- 
served. Note  that  the  region  of  interaction  (inflection  point  ration  in  nega- 
tive slope  side  of  the  power  spectrun)  is  a function  of  the  time  of  interaction 
between  the  different  frequency  components  of  the  packet.  For  a fixed 
receiver  position  as  the  plasma  frequency  is  decreased,  the  packet 
velocity  decreases  thus  increasing  the  interaction  time.  Similarly,  for 
a fixed  plasma  frequency,  the  further  away  from  the  launcher  the  more 
time  the  components  had  to  interact.  This  phenomena  causes  the  in- 
flection point  in  the  power  spectrum  to  move  toward  lower  frequencies. 

The  effects  described  above  were  not  observed  when  the  negative 
pulse  was  used  as  the  exciter.  Even  though  for  this  polarity  Equation 
3-48  predicts  quantitatively  different  results,  the  maximum  amplitude 
available  at  the  plates  and  the  width  of  the  pulse  were  such  that  <r  is  al- 
ways much,  much  less  than  <r  . Figure  4-17  shows  the  plasma  response 

for  the  maximum  pulse  strength  available  with  the  negative  generator. 

Also  shown  is  the  power  spectrum  of  the  response.  As  can  be  seen  from 
the  diagrams,  nonlinear  effects  were  not  observed  using  this  excitation 
as  expected  from  the  values  of  the  nonlinear  parameter  for  this  pulse. 

4.  2 Finite  Axial  Magnetic  Field 


4.  2.  I Linear  Regime 

Since  all  experimental  "idiosyncrasis ''  have  already  been  discussed. 
In  this  and  the  following  subsection  we  need  only  present  the  results. 

Pulse  amplitudes  corresponding  to  field  strengths  of  61  Volts/cm  were 
used  to  excite  the  linear  Bulk  Waves.  The  pertinent  parameters  for  this 
series  of  experiments  were  the  plasma  frequencies  and  axial  field  strength. 
In  figure  4.  1 8,  we  display  time  oscillograms,  recorded  48  cm  away  from 
the  launcher,  as  a function  of  plasma  frequency  with  Bo  fixed  and  in  Fig. 

4.  19  as  a function  of  axial  field  strength  (with  plasma  frequency  fixed). 
Again,  the  region  marked  with  an  A in  these  oscillograms,  cor-es- 
pond  to  electromagnetic  radiation  that  will  still  be  observed  if  the  plasma 
were  not  present.  From  these  diagrams,  we  compute  characteristic 
parameters  that  describe  the  propagating  wave  packet  and  compare  them 
with  the  theoretical  values  obtained  from  aquation  (3.  541. 
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When  the  DC  axial  magnetic  field  is  applied,  the  wavepackets 
are  observed  to  be  noisier  than  for  the  zero  field  case  'compare  the 
oscillograms  in  figure  4-13  to  figure  4-2  for  example).  Therefore 
care  must  be  exercised  in  processing  the  data.  To  calculate  zero 
crossings  and  to  digitalize,  the  waveforms  are  retraced  so  as  to 
smooth  out  the  jitters.  In  fig.  4-20,  the  zero  crossings  of  the  experi- 
mental wavepacket  are  compared  to  those  of  equation  3-54.  For  this 
case,  choosing  the  values  of  and  cu c so  that  one  theoretical  zero- 
crossing is  matched  with  the  corresponding  experimental  value,  all 
other  crossings  are  matched  to  within  10%  of  each  other.  The 
chosen  value  of  up  is  within  20%  of  that  measured  by  the  cavity  and 

w is  within  8%  of  that  measured  using  a Gauss  meter.  Another 
c 

characteristic  of  the  wavepacket  is  the  velocity  of  the  first  peak. 

Fig.  4-21  shows  the  behavior  of  this  velocity  as  a function  of  plasma 
and  cyclotron  frequencies.  The  behavior  is  as  predicted  by  the  theoret- 
ical model  (see  Fig.  3.3).  The  quantitative  depenaance  of  amplitude 
vs',  time  was  not  determined  qualitatively  however,  the  waveforms 
behave  as  expected  from  eq.  3-54. 

For  this  regime,  the  experiments  compared  quite  favorably  to 
the  linear  theory. 

4.  2.  2 Nonlinear  Regime 

As  the  input  electric  field  is  increased  from  61  volts /cm,  gross 
nonlinear  phenomena  begin  to  occur.  From  figure  4-22,  one  observes 
that  signals  arrive  sooner  than  the  linear  waves  and  considerable 
modification  of  the  original  wavepacket  begins  to  take  place.  Close 
examination  of  the  lower  trace,  for  which  the  exciting  pulse  is  largest, 
show  what  appear  to  be  superposition  of  two  wavepackets.  Note  the 
discontinuity  in  phase  between  oscillations.  That  this  indeed  is  so 
can  be  seen  from  the  Fourier  transform  of  the  oscillogram.  The 
spectral  amplitude  and  phase  for  each  waveform  in  figure  4-22  is 
shown  in  figures  4-23  to  4-26.  At  linear  levels  (figure  4-23),  only 
one  peak  (C)  is  shown  corresponding  to  the  wavepacket  of  figure  4- 22a. 
The  slope  of  the  phase  represents  a time  delay  from  some  reference 
point  (the  point  of  arrival  of  EM  signal  for  example).  It  can  be 
used  to  calculate  the  group  velocity  of  the  packet  since  the  distance 
the  signal  has  travelled  is  also  known. 
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At  higher  amplitude  two  distinct  phenomena  occurred  (figures 
4-24  to  4-26).  First  is  the  appearance  of  a low  frequency  pea*c  at 
f = . 14  GHz  (A),  and  second  is  the  splitting  of  the  original  peak  (C) 
into  two  ((B)  and  (C)).  Recalling  the  technique  used  in  obtaining  these 
oscillograms  (loop  gain  constant)  it  is  evident  that  these  changes  are 
due  to  nonlinear  phenomena.  Moreover,  the  possibility  of  these  peaK 
being  higher  order  modes  was  ruled  out  by  the  fact  that  the  group 
delay  for  these  modes  is  much  greater  than  for  the  lowest  order  mode 
(see  equation  3-37). 

From  this  figure  , the  amplitude  of  peak  A is  observed  to  in- 
crease very  rapidly  as  the  input  level  is  increased.  The  frequency 
shift  is  negligible.  No  further  analysis  has  been  made  of  this  packet 
but  it  is  suspected  that  it  belongs  to  the  class  of  solitary  waves. 

The  main  concern  is  with  the  splitting  of  the  packet  (C).  As  the 
input  level  is  increased,  there  is  a shift  to  higher  frequencies  of  the 
original  peak.  A gradual  appearance  of  another  peak  B at  the  low 
frequency  side  indicates  another  wavepacket.  From  the  slopes  or  the 
phase  at  these  two  peaks  the  difference  in  group  velocity  has  been 
calculated.  The  higher  frequency  packet  (C)  was  found  to  nave  a longer 
delay  and  hence  a lower  velocity.  According  to  the  modulation  theory, 
the  difference  between  the  time  of  arrival  is  proportional  to  the 
amplitude  (see  section  3.  3).  This  is  clearly  shown  in  figure  4-27; 
moreover  the  observation  that  the  higher  frequency  wavepacket 
travels  at  a lower  velocity  is  also  in  agreement  with  the  modulation 

theory  results. 

Siijce  there  is  no  absolute  measure  of  the  field  strengths  in  the 
plasma,  a quantitative  comparison  of  this  velocity  difference  between 
theory  and  experiment  is  not  possible.  But  from  the  measured  time 
difference  and  the  group  velocity  of  the  wave,  we  can  calculate  the 
required  field  amplitude.  From  equation  3-63  , to  achieve  the 

largest  observed  separation  of  3.4  ns  over  a distance  of  56  cm,  the 
wave  amplitude  of  the  wavepacket  required  is  about  50  v/cm.  This 
field  represents  a total  energy  of  about  1.  3x10  joules  in  the  wave- 
packet.  The  energy  contained  in  the  exciting  electric  field  impulse  is 
about  3x10*'  joules.  Since  most  of  the  energy  of  the  impulse  goes 
into  the  dummy  load  and  radiated  to  outside  the  parallel  plate  structure 
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it  is  of  interest  to  estimate  the  coupling  efficiency,  i,  e.  how  much 
energy  can  be  coupled  into  the  plasma.  A3  a rough  upper  bound 
estimate,  assume  that  when  the  impulse  field  arrives  at  the  plasma 
tube,  the  2k  V /cm  field  is  applied  to  all  electrons  in  the  region  between 
the  plates  during  the  pulse  duration  (At  - 0.  4 ns),  and  each  electron 
would  gain  a momentum  of  eEAt.  The  total  kinnetic  energy  gained  by 
all  electrons  in  the  region  in  this  way  is  of  the  order  of  9x10”^  joules. 
Comparing  the  orders  of  magnitude  of  the  three  energy  levels,  the 
50  v/cm  field  required  to  produce  the  measured  separation  of  the 
wavepackets  is  certainly  reasonable  for  the  experimental  setup  used. 
The  modulational  splitting  of  the  wavepacket  has  been  observed. 
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V . Concluding  Remarks 


In  the  foregoing,  we  have  described  an  experimental  and  theoretical 
investigation  of  the  propagation  of  finite  amplitude  slow  waves  excited  by 
high  amplitude  baseband  pulses.  The  subject  of  symmetric  slow  mode 
propagation  has  been  incorporated  into  the  framework  of  nonlinear  dis- 
persive waves.  Theoretically,  it  has  been  shown  that  for  typical  experi- 
mental conditions  (i.e.  , = 0,  u^»  w^),  t*le  propagation  of  the 

symmetric  slow  mode  is  described  by  the  Kortewig  de  Vries  equation,  with 
the  magnitude  and  sign  of  the  various  coefficients  of  the  KdeV  equation 
being  different  in  the  various  regimes  considered.  This  implies  that 
given  equal  initial  conditions,  different  effects  in  the  propagation  mode 
will  be  observed.  The  close  interplay  between  theory  and  experiment  leads 
to  a good  understanding  of  the  experimental  results.  Although  linear 
experiments  on  the  propagation  of  this  mode  have  previously  been  done,' 
far  better  resolution  has  been  obtained  here.  To  our  knowledge,  neither 
the  experimental  nor  theoretical  results  of  the  nonlinear  problem  have  been 
previously  reported. 


(31) 


There  are  two  areas  which  need  further  investigation.  The  results 
of  that  investigation  will  further  substantiate  the  theory  developed  here 
and  will  aid  in  the  experimental  verification  of  the  general  theoretical 
results  obtained  for  non-stationary  waves  propagating  in  accordance  with 
the  KdeV  equation 


First,  as  mentioned  in  Section  4.  3.  3,  the  nature  of  the  initial  peak 
in  the  frequency  spectrum  has  not  been  determined  and  second,  nonlinear 
experiments  where  the  polarity  of  the  exciting  pulse  is  negative  must  still 
be  performed.  We  will  not  consider  the  first  area  further  and  instead 
comment  only  on  the  second  problem.  The  first  implication  of  reversing 
the  polarity  is  that,  theoretically,  the  phenomena  observed  for  the  weak 
field  case  should  now  occur  for  the  zero  magnetic  field  case  and  vice  versa. 
This  is  a consequence  of  the  sign  of  the  nonlinear  term  in  equations  (3-47) 
and  (3-48). 

In  conclusion,  it  must  be  noted  that  in  order  to  experimentally 
obtain  a wider  range  of  parameters,  the  experiment  as  described  should 
be  modified  in  two  ways:  the  duration  of  the  pulse  should  be  increased, 
and  the  cathode  arrangement  changed.  Increasing  the  pulse  duration 
enhances  the  nonlinear  effects  (see  equations(3 -57)  and  (3-58)),  whereas 
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changing  the  cathode,  to  a flat  cathode  for  instance,  will  eliminate  the 
rotational  instability  of  the  column,  observed  at  high  magnetic  fields. 
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VI.  Appendices 

6.  1 Formal  derivation  of  operators 

Equation  (3-10)  may  be  written  approximately  up  to  second 
order  in  the  field  variables  as: 


u - 


(6-1) 


where, 

md.u  ^ + mi/  u^3-enE-euu^xH  (6-2) 

t—  c - o — - -o 

and 


mdfu<2)  + mi/u<2>  - .a^.u^u^-wE 
t~  c—  n — — — 


Formally  solving  (6-2)  and  (6-3), 


(1) 


-1 


(md.  1 - eu  H x 1 + mi/  1)  en  E 

~Q  ~ O S Css  O — 


C* 


eu  u(2)  x H 

o-  -o 


(6-3) 

(6-4) 


u(2)  * - (mdl  - eu  H xl+mi/1)  1 (r^7.  u^V'^+enE]  (6-5) 

t ss  0 — 0 a Cs  L n — — — 

O 

Defining  the  inverse  operator: 

n(dt'90'l/c)3  (mdti  *efio90xi't'mi/ci  )_leno  (6'6) 

the  total  current  from  (6-1)=  (6-4)-(6-6)  aid  (3-13)  is  given  by: 

u 3 - n { E + -TT-  V.  it  E tr  E - — E V.  E } (6-7) 

— ~ _&  a —a  — en  — — 

n e o 

o 

which  is  used  to  arrive  at  equations  (3-1 5). 

Defining  the  time  Fourier  transform  of  a field  variable  as: 


f(r.t) 


r..  . -iwt  dw 

-0C 


Applying  the  above  to  equations  (6-6)  and  (6-7),  we  obtain: 
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II  + ~ ~ Qt'J,  H ,v  ) p < b 

- € -J  - O C 

O 

(6-11) 

I P > b 

and  the  longitudinal  field  components  may  be  obtained  from  the  trans- 
verse components  via: 


Note  that  since  the  operator  Lq  takes  different  forms  in  the  range 
0 < p < oo,  likewise  £a's  will  have  different  representations. 

In  general,  the  total  field  corresponding  to  and  H will 

have  both  Eff  and  components.  This  impJies  that  no 

TE  or  TM  mode  representation  is  possible;  therefore  a complete  set 

of  modes  must  be  composed  of  TE  and  TM  modes.  In  the  following 

analysis,  we  will  restrict  to  the  cases:  a)  zero  static  magnetic  field, 

b)  infinitely  large  magnetic  field  and  c)  weak  magnetic  fields,  i.  e. 

(*>te<u>pe,  and  the  column  is  surrounded  by  a perfect  condition.  The 

first  two  cases  result  in  substantial  simplification  in  the  form  of  the 

eigenvectors  and  eigenvalues.  For  these  cases,  Eqs.  (2-3)  can  be 

satisfied  with  either  E or  H equal  to  zero.  This  implies  that 

az  az 

two  independent  solutions  are  possible.  Since  neither  the  nonlinearity 
nor  the  boundary  couple  the  two  solutions  and  since  initially  only  E 
modes  are  generated  by  the  electric  current  source,  only  these  will 
be  considered  in  the  field  representation.  For  the  low  magnetic  field 
case,  we  make  the  assumption  that  the  waves  are  slow  and,  therefore, 
are  primarily  TM  type,  i.  e.  , the  actual  field  eigenvectors  are  approx- 
imated by  quasistatic  eigenvectors. 

An  implicit  assumption,  when  representing  the  total  fields  as  a 
superposition  of  modal  fields,  is  that  the  modal  set  (or  eigenvectors) 
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is  complete,  i.  e.  , they  form  a basis  in  terms  of  which  any  arbitrary 
vector  may  be  represented.  A statement  of  completeness  may  be 
expressed  as  follows 


I K* 

at  t 


(O1) 


i s(n  - q«  ) 


(6-13) 


The  sumation  in  (6-13)  is  to  be  interpreted  as  an  integral  for  the  case 
of  continuous  spectrum  and  as  a sum  for  the  discrete  spectrum.  A 
direct  procedure  for  determining  the  spectrum  utilizes  the  character- 
istic Green's  function  for  E modes.  The  properties  of  the 

Green's  function  in  the  complex  k plane  are  investigated;  and  from 
pole  and  branch  singularities,  which  contribute  to  the  integral  of 
Green's  function  along  a suitable  contour  in  the  complex  k plane, 
the  eigenvectors  are  obtained.  In  this  case,  the  question  of  complete- 
ness is  assured  although  normalization  as  in  (3- 19a)  is  not.  This 
procedure  will  be  used  subsequently  for  the  determination  of  the  mode 
vectors  of  the  zero  axial  magnetic  field  case. 

For  the  other  two  cases,  we  proceed  directly  to  obtain  that  part 
of  the  spectrum  associated  with  guided  waves. 


6.  2.  1 Eigenvectors  for  zero  magnetic  field 

When  H is  zero,  then  tensor  e degenerates  into  a scalar 
quantity.  In  the  plasma  region,  e — • 1(1-  — Hr- ) . Since  for  an  E 
mode  . zqx  E^  = 0,  we  have  (from  equation  (6-10)  for  the  scalar 
Green's  function:  (40) 


T • vtg  + eg  = 


S(p-p') 


(6-14) 


where 


2 2 
U>  U«  - € K 

no  o a 

2 _ 2 
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p < b 
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Boundary  conditions  for  g may  be  obtained  from  (6-14),  since  it  is 
valid  for  all  p.  By  integrating  across  p = b and  p = p',  we  have 
the  conditions: 


1)  g is  continuous  at  p = b and  p = p 1 

2)  Vtg  /°*  - - 1 

K ' ' 


(6-15) 


3)  V /bb  * 0 


for  a symmetric  source  of  electric  current  located  at  p 
solution  to  (6  - 1 4)  satisfying  (6 - 1 5 ) is: 


= b , the 


■J  (K  b)  J (K  p) 
° M ° X1 


* 


-JJk  b)  H (K  b)  H (/C  P) 

o i2  ° A2 


° x 


1 


H“(k  b) 


& 


x2 


o < b 


p > b 


(6-16) 


€ J b)  (k  b)  J_(k  b)  H,U  b) 


where  a = 


1 


° X 


A1 


(6-17) 


and  J.  , J , H , H.  are  Bessel  functions  of  the  first  and  second  kind, 

1 Or  O L 

respectively. 

The  spectral  representation  in  (6-13)  may  now  be  obtained  by 
integrating  the  scalar  Greens  function  along  a contour  that  encloses  all 
singularities  of  g.  The  desired  relationship  is  given  by^^: 
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t2 
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(6-13) 


U b)  H (K  b)  H U P)d  k 


o x 


the  transverse  component  of  E is  obtained  from  <p  using  , 


E = 
*at 


2 

k ue 

•4.  0 


7t8 


i => 


p < b 
p > b 


(6-19) 


The  integrand  in  (6-18)  is  analytic  everywhere  except  for  a finite 

number  of  poles  and  branch  points.  For  real  < , only  a single  pole 

exists  at  & = 0;  the  branch  point  corresponds  to  k =0  for  which 

~ -1-2 

<_  s + m V . The  pole  contribution  gives  rise  to  a guided  surface 


a - 


o o 


mode,  while  the  branch  cut  contribution  gives  rise  to  a radiation  field. 

For  the  excitation  problem  both  contributions  must  be  included, 

specially  at  the  source.  In  the  far  field  we  only  consider  the  guided 

wave,  but  as  observed  in  Chapter  III;  we  do  not  know  the  absolute  field 

strength  of  the  guided  mode  since  we  negiect  to  include  the  amount  of 

energy  extracted  from  the  source  by  the  other  mode  types.  For  guided 

modes,  k = i[/<  I 

4 ' -4 

Lode 


(46) 


guided  modes: 


< = i | k | and  from  (6-18),  we  have  for  the 
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with  given  by: 
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1 xi  ° x2 
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From  (6-19)  the  transverse  field  components  are  given  by: 
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(6-22) 


Even  though  the  potential  functions  are  orthonormal:  the  transverse 
components  of  the  eigenvector  do  not  satisfy  condition  (3-19a).  Proper 
normalization  is  introduced  by  defining  a normalization  factor  as  follows 

e£(x  p) 

1 1Z 

———  -P  dp 

x *"wt  Jd  A^ 
x2  oV  Ka 


the  above  can  be  rewritten  as: 
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where 
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(6-24)  is  evaluated  explicitly  by  using  the  formulae 


fl  U p)pdp  = -i  b2{l*  U b)-(l+  -;-1  x )lhK  b)} 
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the  characteristic  equation  (6-20),  and  the  approximate  disr.  : - ion 

equation  (3-15).  In  the  region  < b < 1 , we  can  expand  the  fsessel 

1 1 

functions  in  the  small  argument  regime  and  obtain: 


P 

where 


2 2 
u v , 
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and  C = speed  of  light. 


with  vq  given  by  (3-34) 


The  longitudinal  component  of  the  eigenvector  as  defined  in  (3-22)  is 
obtained  from  (6-23),  (6-21)  and  (6-12). 

6.  2.  2 Eigenvectors  for  large  Magnetic  fields 

In  the  limit  of  large  magnetization,  i.  e.  , u /<*)  » 1 u »u  ; 

Cm  C P 

transverse  electron  motion  is  effectively  inhibited  and  the  plasma  be- 
comes non-gyrotropic.  The  tensor  e in  equation  (6-10)  for  p < b 


- LOO- 


The  guided  eigenvector  components  may  be  obtained  from  a 
potential  $ as: 

5.  = - —z  * (6 


H = z x E 

- a.  k -o  -a, 
t a t 


where  <b  obeys  the  equation: 
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equation  (6-27)  is  solved  subject  to  two  different  boundary  conditions: 
1)  column  surrounded  by  a perfect  conductor  and  2)  column  in  free 
space. 

1)  for  a perfect  conducting  boundary  the  solution  with  proper 
behavior  at  p = 0 and  p = b is  : 

<t>~  A J (k  p)  (6-30) 


with  k b = p > where  p is  the  n zero  of  the  Bessel  function  J . 
x ^ rn  rn  o 

Using  (6-26)  and  (6-12),  the  eigenvector  components  are  given  by: 
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the  normalization  constant  is  determined  from  the  orthonormality 
condition  (3 -19a): 
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real  , we  have 


If  the  a index  is  separated  into  two  indecis:  m = mode  type  , 
n = mode  in  mode  type;  we  have: 

d)  — il;  —l;  since  only  E modes  are  considered. 
a mn  n 

2)  Column  in  free  space 

The  solution  (6-30)  of  (6-28)  is  still  valid  in  the  region 
p < b,  the  plasma  region.  For  p > b,  the  solution  that  satisfies 
boundary  conditions  at  p = b and  appropriate  to  the  discrete  spectrum 
is: 

J U b> 
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where  k and  < are  real  and  given  by  (see  6-29)  : 
x 2 X1 
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and  the  boundary  characteristic  equation 
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from  (6-30),  (6-31)  and  (6-26);  the  field  components  of  the  eigenvector 


are: 
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using  above,  we  can  expand  the  collision  coefficient  as  in  (3-21) 
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The  integrals  in  the  bracket  are  evaluated  using  the  formulae 
in  section  6.  2.  After  approximating  in  k b < 1 region: 

2 2 

• .2  „ tx  I“(k  b) 


XU 


= -f—r  2#- 

(u  +VT) 


■1 


using  the  expression  for  (3  from  6.  2 and  the  linear  approximation  for 
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independent  of  frequency. 


(6-35) 
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for  p < b. 
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for  p > b,  solution  (6-31)  holds.  The  normalization  constant  A'  is 
obtained  via  (3- 19a)  as  follows: 
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evaluating  the  integral  as  in  (6-24);  the  approximate  value  for  the 
normalization  constant  is: 
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(6-34) 


the  desired  eigenvector  is  given  by  (6-31),  (6-3  3)  and  (6-34). 

6.  3 Explicit  evaluation  of  nonlinear  and  collision  operators 

Using  the  eigenvectors  derived  in  the  previous  section  and  the 
approximate  dispersion  relations  presented  in  section  3.  2,  we  can 
explicitly  evaluate  the  collision  and  nonlinear  operators  in  equation 
(3-28). 

6.  3.  1 Zero  Magnetic  field 

When  Hq  = 0,  the  operator  £ in  equation  (3-21)  takes  the  form: 
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using  above,  we  can  expand  the  collision  coefficient  as  in  (3-21) 


V(y)  . (+a.v+<|)  .1  / 5a-  E*  pdp 


Id  + , 


. VV‘  d k'oV 

,,2,2  u2fi  U 2 pdp  +i  5 


u + v 


pdp  } 


A1 


The  integrals  in  the  bracket  are  evaluated  using  the  formulae 
in  section  6.  2.  After  approximating  in  < b < 1 region: 
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the  above  equation  is  the  same  equation  obtained  by  Trivelpiece  et  al/ 

9 K 

since  — - — - v for  u « u 
o p 

Neglecting  the  effects  of  collisions,  the  nonlinear  term  in  equation 
(2-21)  is  expanded  as: 


i[Efl(w).N'(«aEa(u))] 


^a(‘-“l)^(*-u2)V  ??(?aL  • 5p“ 

+ 4a(z,ul)dz4a(z.'J2)EaZ(Ea1  . Epu)] 


u, 


2 

eu  e 


«i) 


4 -2-2-  C(  (z.u,)  [E  • (z  3 tf(2.u,)E„ 
u.u.m  a 2 L -a  ~o  z tj  1 -a 


1 2 


+ (*.«!>  • VtE^L]  . E*“}  6(u-u1-u2)du1du2 


ui  . 


where  Ea  is  the  eigenvector  associated  with  the  amplitude  m^(z  , un).  In 
order  to  simplify  the  above  equation,  we  expand  CC  (z,u^)d^Cf^(z,u^)  up  to 
second  order  as: 


tfa(z,ul)8z^(*,u2)  =^a(z>u1)  Ua(u2)*a(z,u2) 


(6-36) 


Using  this  expansion,  the  first  square  bracket  term  becomes  zero 
while  the  second  is  expanded  to: 


ee  w2  „ CL  (z,u,  )C(  (z,u,)  b u,  u.  u*  o u, 

-2-2-  f-S. - — 3 Mj[E  29  E 1 E 4E29E1E 

mu  J u^  u^  £ ap  p ap  ap  ap  p ap 


az 


U2  W.  U*  fc>-  w.  u* 

4 i k (u. )E  ^ E 1 E 4i  K (u,  )E  “,E  1 E ] pdp}6(u-w,  -uOdu-du, 
a 1 az  ap  ap  a 1 az  az  azJ  r 12  12 


(6-37) 


putting  in  for  the  eigenvector  components  from  section  6-2  we  obtain 
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equation  (3-29). 

Putting  in  for  (3  from  section  6.  2 and  letting  I.  (*  p )-  -s  < p » 

1 - 1 x 1 

(2)'  (2)  / 1 
I,  -<  /2  I (it  p ) i.  e.  < b<l  regime,  the  nonlinear  term 

1 M 0 J-l  X1 

corresponding  to  P is  evaluated  as: 


u.uyj  b 

= m f-TTZ 3TTrta(i,Ul)  (*.u2)6(u-«1-w2)duldu2  (6-38) 

2 (u  v,v  be 
v p d o o 


where 


vd  = 


(1-v^/c2) 


and  for  <a(u)  the  linear  contribution  has  been  used.  Similarly  the  con- 
tribution due  to  u is  evaluated 


«ff  “i  “ 

" ‘ mJ  1 . .6  6 6~§  T (1).  (2).  (o)  l (1) 

(ebtjvv,)4  I'  'I'  'I  ' 'o  k k 
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P dp 
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pod' 


•1  X1 


V Id- 

o 2 


, „6  6 4 ! (Dr  (o) 

(v  b v , e ) o o o 
o do 


b j(1)j(2)j(o) 


A1 


^ (z.WjJ^z.UgJSfu-Uj-w^  dw1du2  (6-39) 

K 

1 1 

letting  1^  ~ — ^ — Iq  and  expanding  Iq  to  first  order  (for  k b < 1 ). 

the  integrals  in  (6-39)  can  be  readily  evaluated.  Combining  the  result 
with  equation  (6-38),  we  arrive  at  equation  (3-39). 

6.  3.  2 Infinite  Magnetic  field. 

When  the  magnetic  field  is  large,  the  nonlinear  term  in  equation 
( 3-21)  reduces  to: 
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+ ]fEa  ^a^a  pdp 

14  o z z z 


b u2  oj.  u* 

4 (z,u.)  4?  (z,u,)  f E “Vt'  E-  1 E pdp}  6((j-u)..u,)du).dij,  (6-40) 

a i a c j a i — a—,  a i c i l 


T z 


from  equation  (3-13)  and  (3-10),  we  have  to  first  order: 

2 

T u z 


Using  above  equation  and  the  expansion  of  (6-36)  in  (6-39),  we  obtain 
for  nonlinear  term: 
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Putting  in  for  E from  section  6.2, 
z 


= - ^rf{(u :5;  + “I-)  Ka(“1)^.(2.u1)^(z.U2,A,(0J)A'(0J  )A'(U  ) 
1 2 


f Jd  (<  p )pdp  } 5 (cj-u. -oj-)doj,du0 
J0  o Xl 


(6-41) 


the  form  of  (6-41)  is  the  same  for  both  cases:  column  in  free  space  and 
column  surrounded  by  perfect  conductor.  The  difference  lies  in  the 
explicit  form  the  normalization  coefficient  A' (u)  takes.  When  the  column 
is  surrounded  by  a perfect  conductor  the  normalization  is  (from  section 
6.  2):  k 

A'(u)=  (6-42) 


(u 


^ M 

KaO*U»  )bJ,U 

a 0 u2  1 M 


b) 


where  k b = p . p is  the  ntd  zero  of  J . Assuming  that  only  the 
j n n o 
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lowest  order  mode  is  present  (higher  order  modes  are  slower  and 
consequently'  attenuated  at  a faster  rate)  3"^, 


< b = 2.  40  5 
A1 


and  equation  (6-41)  becomes: 


(2u,+u,)  u ,,  Anc\l  -L  ~o 


'1  2;  , _x_  2 (2.  40 5)  _o 


/ P )pdp 


1 


~4  - —5  T 3, 

u)_  o b J,  k b 

P 1 XI 

^*.U2)  Mu-id^  -id^dw^du^ 


} (z  , W ) 
<2  i 


(6-43) 


In  the  above  derivation,  it  is  assumed  that  id  < <d  , so  that  the  linear 

P 

portion  of  the  dispersion  equation  gives  the  maximum  contribution  to 
the  nonlinear  term.  For  the  column  in  free  space, 


A'(u)  = 


u bJ  (k  b) 


p o'  x 
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(6-44) 


and  (6-41)  becomes 
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(6-45) 


where  the  same  assumptions  leading  to  (6-43)  have  been  used.  (6-43) 
and  (6-45)  lead  to  (3-42). 

The  collision  operator  for  large  magnetic  fields  is  given  by: 


v = L - L = i 
= o 


where 
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The  collision  term  takes  the  form  (for  v / y < 1), 


('-•  . X 40  = i J E 

<2  — d J 


* u;  v e 

a V <“l^l>PdP 

2 2 U rvc 


putting  in  for  the  longitudinal  component  of  the  eigenvector, 


= i*[  Ja* 2('j)  p)pdp] 


UDVce0 

~2~~l 

(J  + v 


(6-46) 


Again,  the  above  form  applies  for  both  column  in  free  space  and  column 
surrounded  by  conductor.  Using  (6-42),  the  collision  coefficient  when 
the  column  is  surrounded  by  a perfect  conductor  is. 


y(u)  = 


for  v /<j  < 1 


and  again  we  have  assumed  that  u < , so  that  the  major  contribution 

to  k arises  from  the  linear  term. 
a 

Similarly,  for  the  column  in  free  space: 


Y(u)  = 


,2  T 2t 


X1 


b) 


)pdp 


6.  3.  3 Weak  Magnetic  field. 

To  obtain  the  Kernel  given  in  (3-30);  the  nonlinear  term  in 
equation  (2-21)  is  expanded  in  component  form,  and  consistent  with  the 
slow  mode  assumption,  we  only  retain  the  term  arising  from  the  coupling 
of  the  longitudinal  electric  field.  For  (j  < uc  < , this  term  is  given 

by: 
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Using  the  expression  for  the  z component  of  the  eigenvector  (from  6.  2.  3) 
to  explicitly  write  the  inner  product: 


2 2 

(J  w 
D C 


— - f[  (2 ? r ) 2 — + 5 — 3 r~]  <■  (u.K?  (z.w.)Z?  (z.u,) 

m v 2.  2 u.U-U  2,  2,  2 ala  i a £ 
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pc  1 p c 

A(uL)  A(u2)A(u)  fj*  (K^^p)pdp  5 (u-u1-u2)du1du2 


Putting  in  for  the  normalization  coefficient  A(u)  and  proceeding  as  in 
6.  3.  2 . we  arrive  at  the  final  form  (3-44). 
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